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PREFACE 


This  report  was  prepared  by  Thaddeus  C.  Johnson,  Research  Civil  Engineer,  Civil 
Engineering  Research  Branch;  Richard  L.  Berg,  Research  Civil  Engineer;  Edwin  L. 
Chamberlain  and  David  W.  Cole,  Research  Civil  Engineers,  Applied  Research 
Branch,  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory. 

The  study  was  funded  under  DA  Project  4A762730AT42,  Design,  Construction 
and  Operations  Technology  for  Cold  Regions  under  several  work  units.  The  Federal 
Highway  Administration  (FHWA)  and  the  Federal  Aviation  Administration  (FAA) 
shared  equally  in  the  funding,  under  FHWA  Order  No.  8-3-0187,  Full-Scale  Field 
Tests  to  Evaluate  Frost  Action  Predictive  Techniques. 

The  research  findings  summarized  here  are  the  product  of  a  large-scale  effort  that 
required  the  participation  of  several  agencies  and  many  persons.  Particularly  signifi¬ 
cant  were  the  research  contributions  of  Dr.  Gary  Guymon,  University  of  California 
at  Irvine,  who  played  a  key  role  in  developing  the  frost  heave  model,  and  Dr.  Lynne 
Irwin,  Cornell  University,  who  developed  the  pavement  response  model  NELAPAV 
and  made  other  contributions  to  the  research  on  thaw  weakening.  The  following 
were  members  of  a  panel  of  consultants  who  provide  guidance  and  periodic  review 
of  the  work:  Dr.  B.  J.  Dempsey,  University  of  Illinois  (1978-79);  Dr.  D.G.  Fredlund, 
University  of  Saskatchewan;  Dr.  M.E.  Harr,  Purdue  University;  Dr.  R.D.  Miller, 
Cornell  University  (1981-84);  Mr.  E.  Penner,  National  Research  Council  of  Can¬ 
ada;  Dr.  M.W.  Witczak,  University  of  Maryland. 

Many  persons  at  CRREL  participated  in  the  research  effort.  Particularly  note¬ 
worthy  are  the  technical  contributions  of  D.  Bentley,  D.  Carbee,  G.  Durell,  J.  Inger- 
soll,  D.  Keller,  J.  Major  and  R.  Roberts. 

The  authors  greatly  appreciated  the  permission  of  the  Massachusetts  Department 
of  Public  Works  and  the  Albany  County  Airport  to  use  paved  facilities  under  their 
jurisdiction  as  test  sites.  Mr.  Dennis  Rice,  then  of  the  MDPW,  and  Mr.  Jack  Masko, 
Airport  Manager,  were  especially  helpful. 

And  finally,  the  authors  enjoyed  strong  and  unfailing  support  from  the  sponsors 
throughout  this  work.  Principal  technical  liaison  was  provided  by  Mr.  A.  DiMillio 
of  the  Federal  Highway  Administration,  and  the  authors  particularly  appreciate  his 
active  participation.  Mr.  H.  Tomita  of  the  Federal  Aviation  Administration  was 
consistently  helpful,  and  Mr.  A.  Muller,  representing  the  Corps  of  Engineers,  ac¬ 
tively  supported  the  research  and  its  sponsorship  by  the  Corps  from  its  inception. 

This  report  was  technically  reviewed  by  F.H.  Sayles  of  CRREL. 
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Frost  Action  Predictive  Techniques  for  Roads  and  Airfields 

A  Comprehensive  Survey  of  Research  Findings 


T.C.  JOHNSON,  R.L..  BERG,  E.J.  CHAMBERLAIN  AND  D  M.  COLE 


INTRODUCTION 

Six  years  of  intensive  research  has  significantly 
advanced  the  state  of  knowledge  and  the  ability  to 
predict  the  effects  of  frost  action  on  pavement 
performance.  The  two  principal  adverse  effects  of 
frost  are  ice  segregation,  causing  heave  and  tran¬ 
sient  pavement  roughness,  and  thaw  weakening  of 
subgrade  and  unbound  base  materials,  causing  ac¬ 
celerated  pavement  cracking  and  pavement  defor¬ 
mation  (Johnson  et  al.  1975).  In  the  past,  design¬ 
ers  of  new  and  rehabilitated  pavement  projects 
have  necessarily  applied  thickness  design  methods 
that  are  essentially  empirical,  reflecting  pavement 
substructure  designs  that  have  performed  ade¬ 
quately  under  similar  conditions  in  the  past.  While 
the  empirical  approaches  have  been  quite  success¬ 
ful,  it  has  long  been  an  important  goal  not  only  to 
improve  them  but  specifically  to  develop  quantita¬ 
tive  methods  for  predicting  the  surface  heave  of  a 
pavement  section  and  for  evaluating  the  seasonal 
changes  in  supporting  capacity  of  subgrade  and 
base  materials  that  would  affect  pavement  perfor¬ 
mance  under  traffic  loads  (Johnson  et  al.  1975). 

Having  a  common  interest  in  these  goals,  three 
agencies  jointly  sponsored  the  research:  the  Fed¬ 
eral  Highway  Administration,  the  Federal  Avia¬ 
tion  Administration  and  the  U.S.  Army  Corps  of 
Engineers  (DiMillio  and  Fohs  1980).  The  research, 
which  spanned  the  years  from  late  1978  through 
1984,  included  equipment  development,  field  and 
laboratory  experiments,  and  development  of 
mathematical  models.  The  investigations  focused 
on  four  principal  study  areas: 

•  Selection  and  validation  of  the  most  effec¬ 
tive  laboratory  index  tests  to  indicate  the 
susceptibility  of  soils  to  detrimental  frost 
action. 

•  Development  of  a  soil  column  device  with 
provisions  for  nondestructive  monitoring  of 
changes  in  soil  moisture  content  and  density 
during  freezing  and  thawing. 


•  Improvement  and  validation  of  a  mathe¬ 
matical  model  of  frost  heave  that  had  been 
developed  earlier,  and  incorporation  of 
processes  that  take  place  during  and  after 
thawing. 

•  Development  of  laboratory  test  methods 
for  characterizing  seasonal  changes  in  the 
resilient  modulus  of  a  wide  variety  of  gran¬ 
ular  soils,  and  validation  of  these  methods 
by  means  of  in-situ  deflection  testing  of 
pavements. 

The  research  also  included  consideration  of  means 
of  implementing  the  findings  in  pavement  design 
practice,  and  preparation  of  flow  charts  showing 
how  this  might  be  accomplished. 

The  research  was  performed  by  the  Army  Cold 
Regions  Research  and  Engineering  Laboratory 
(CRREL).  Key  investigations  were  conducted  by 
consultants  on  certain  phases  of  the  work,  and  we 
also  assembled  a  board  of  general  consultants  who 
provided  guidance  and  periodic  review  of  the  ac¬ 
complished  work.  The  laboratory  testing  and 
analysis  were  conducted  at  CRREL’s  facilities  in 
Hanover,  New  Hampshire,  and  field  sites  at  Win- 
chendon,  Massachusetts,  and  Albany  County  Air¬ 
port,  New  York,  were  used  for  in-situ  testing  and 
data  collection. 

This  report  presents  a  comprehensive  survey  of 
the  research  findings  in  each  of  the  four  study 
areas.  Details  of  the  investigations  have  been  doc¬ 
umented  in  one  or  more  research  reports  on  each 
topic,  for  which  citations  are  given  throughout 
this  report. 

FIELD  TEST  SITES 

Field  test  sites  were  needed  to  serve  as  a  source 
of  samples  of  subgrade  soils  and  base  materials 
for  roads  and  airfield  pavements,  and  as  test  beds 
where  the  performance  of  these  soils  and  materials 
in  pavement  sections  could  be  monitoreu  and  test- 
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Figure  I .  Winchendon  test  site. 


ed  under  varied  conditions  of  temperature,  mois¬ 
ture  and  freeze-thaw  action  and  under  applied 
loads.  Two  sites  were  used  for  these  purposes. 

The  first  is  the  Winchendon,  Massachusetts, 
test  site,  constructed  in  1978  by  the  Massachusetts 
Department  of  Public  Works  (MDPW).  The  site 
consists  of  24  soil  test  sections  (Fig.  1).  Twelve 
soils  are  used,  each  in  embankments  of  two 
heights.  The  higher  embankments  were  chosen  for 
the  six  test  sections  used  in  this  research  The  sec¬ 
tions  consist  of  about  50-90  mm  of  asphalt  con¬ 
crete  and  1.5  r.i  ot  test  soil  (either  Ikalanian  sand. 
Graves  sand.  Hart  Bi others  sand,  Hvannis  sand, 
dense-graded  stone  or  Sibley  till)  overlying  the 
natural  subgrade,  a  dean  gravelly  sand  (see  grain- 
size  curves,  Fig.  2).  The  water  table  is  about  1 .4  m 


below'  the  pavement  surface.  Data  were  collected 
from  each  of  the  six  test  sections  (Johnson  et  al. 
1986a,  Guymon  et  al.,  in  prep.)  during  the  freeze- 
thaw-recovery  seasons  of  1978-79  and  1979-80, 
including  the  following: 

•  Temperature,  monitored  by  thermocouples 
at  the  surface,  within  the  asphalt  concrete, 
and  at  various  depths  in  the  subgrade  (Fig. 

3). 

•  Depth  of  frost,  monitored  by  thermo¬ 
couples  and  by  frost  tubes  previously  in¬ 
stalled  by  MDPW  (Fig.  4). 

•  Depth  to  the  water  table,  monitored  in  ob¬ 
servation  wells  previously  installed  by 
MDPW. 

•  State  of  stress  in  soil  moisture,  monitored 
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Figure  2.  Grain-size  distributions  of  Winchendon  test  soils  and  natu¬ 
ral  subgrade. 


Figure  3.  Ground  temperatures  prevailing  during  plate-bearing  tests  in  the 
Ikatanian  sand  test  section. 


by  means  of  moisture  tensiometers  installed 
in  the  subgrade  at  depths  of  152,  305,  610 
and  914  mm  below  the  paved  surface  (Fig. 
4). 

«  Vertical  displacement  of  the  paved  surface 
under  load,  measured  on  the  paved  surface 
at  a  central  point  on  each  of  the  six  sections 
on  13  occasions  by  repeated-load  plate¬ 
bearing  (RPB)  tests  and  by  falling-weight 
deflectometer  (FWD)  tests  on  nine  occa¬ 
sions  (Fig.  5). 

•  Vertical  displacement  of  the  paved  surface 
caused  by  freezing  and  thawing,  measured 


by  level  surveys  conducted  by  MDPW  (Fig. 

4). 

•  Meteorological  data,  including  precipita¬ 
tion  and  air  temperature,  monitored  by 
means  of  a  precipitation  gauge  and  a  ther¬ 
mograph. 

The  second  test  site  is  at  Albany  County  Air¬ 
port,  located  .0  km  north  of  Albany,  New  York, 
where  two  pavements  were  selected  for  field  test¬ 
ing:  Taxiways  A  and  B  (Fig.  6).  Taxiway  A  is  a 
new  pavement  that  was  under  construction  in  1979 
when  we  installed  the  instrumentation,  and  the 
first  loading  tests  were  performed  prior  to  any  ait 
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Figure  4.  Frost  depth,  heave  and  moisture  tension  in  the  Graves  sand 


test  section. 


craft  traffic.  The  pavement  cross  section  (Fig.  7) 
consists  of  330  mm  of  asphalt  concrete,  584  mm  of 
crushed  stone  base  and  914  mm  of  gravelly  sand 
subbase  over  a  subgrade  of  silty  fine  sand.  The 
grain-size  distribution  curves  (Fig.  8)  show  that 
both  the  base  and  subbase  materials  have  signifi¬ 
cant  fine  fractions,  with  about  12°/o  passing  the 
No.  200  sieve  and  about  997o  finer  than  0.02  mm. 
In  both  cases  the  fines  are  nonplastic. 

Taxiway  B,  constructed  many  years  ago  (possi¬ 
bly  in  the  1940s),  is  no  longer  used  as  a  taxiway 
and  is  used  only  occasionally  for  parking  light  air¬ 
craft.  The  surface  is  uneven  and  the  asphalt  con¬ 
crete  is  aged  and  severely  cracked.  The  asphalt 
concrete  is  76  mm  thick  and  overlies  about  102 
mm  of  asphalt-penetration  macadam  stone  base 
and  127  mm  of  gravel  subbase  (Fig.  7).  The  sub¬ 
grade  is  silty  fine  sand.  The  penetration  macadam 
base  is  aged  and  deteriorated  to  the  point  of  being 
quite  friable  and  is  not  clearly  distinguishable 
from  the  subbase.  Consequently  the  samples  of 
the  subbase  whose  composite  grain-size  distribu¬ 
tion  is  shown  in  Figure  8  were  taken  to  represent  a 
single  229-mm  base-subbase  layer. 

Boring  records  showed  groundwater  about  2.3 
m  below  the  pavement  of  Taxiway  A  and  1.8  m  at 
Taxiway  B.  We  installed  groundwater  observation 
wells  at  both  taxiways  and  found  that  the  depth 


fluctuates  seasonally  between  about  1 .5  and  2.0  m 
at  both  sites. 

The  following  data  were  collected  at  each  of  the 
test  pavements  at  Albany  County  Airport  (Guy- 
mon  et  al . ,  in  prep.,  Johnson  et  al.  1986b)  during 
the  four  freeze-thaw-recovery  seasons  from  1979 
to  1983,  except  as  noted: 

•  Temperature,  monitored  by  thermistors  at 
the  surface,  within  the  asphalt  concrete, 
and  at  various  depths  in  the  base,  subbase 
and  subgrade  (Fig.  10). 

•  Depth  to  the  water  table,  monitored  in  ob¬ 
servation  wells  adjacent  to  the  taxiways. 

•  State  of  stress  in  soil  moisture,  monitored 
by  tensiometers  at  depths  of  about  533,  813, 
1118  and  1 524  mm  below  the  pavement  sur¬ 
face  at  Taxiway  A  and  480,  640,  840  and 
1140  mm  below  the  pavement  surface  al 
Taxiway  B  (Fig.  1 1 ). 

•  Vertical  displacement  of  the  pavement  sur¬ 
face  under  load,  measured  by  the  RPB  de¬ 
vice  on  two  occasions  and  by  the  FWD  de¬ 
vice  on  10  occasions  in  the  1979-80  season, 
and  1 1  times  by  the  FWD  device  tn  1982-83 
(Fig.  12). 

•  Vertical  displacement  of  the  pavement  sur¬ 
face  due  to  freezing  and  thawing,  measured 
bv  rod  and  laser  level  survevs. 
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Figure  6.  Test  site  locations ,  Albany  County  Airport ,  Albany,  New  York. 
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•  Meteorological  data,  including  precipita¬ 
tion,  wind  speed  and  direction,  and  air  tem¬ 
perature,  obtained  from  the  National  Wea¬ 
ther  Service  station  at  the  airport. 

In  addition  to  the  data  collection  and  in-sit u 
testing  at  the  Winchendon  and  Albany  County 
Airport  sites,  we  also  used  samples  from  the  same 
test  sections  for  related  laboratory  tests.  Core 
samples  of  the  asphalt  concrete,  102  mm  in  diam¬ 
eter.  were  taken,  and  5'  mm-diameter  undis¬ 
turbed  samples  ot  the  finer  soils  were  obtained  in 
the  tall  prior  to  tree/e  up  Once  frost  had  ad¬ 
vanced  to  MiHisient  depth  m  the  sections,  core 
sample'  ot  the  tro/en  soil  were  taken,  except  those 
materials  containing  numerous  graxel-si/e  trag 
merits  Hulk  samples  ot  about  40  M)  kg  also  were 
obtained  t'oin  ea  i  -oh  mJ  base  material 
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Figure  7.  Pavement  sections.  Taxi  wavs  I  ami  H 
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Figure  8.  Grain -size  distribution  of  materials  from  Taxi  ways  A  and 
B. 
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Figure  10.  Ground  temperatures  during  plate-bearing  tests  at  Taxiway  A, 
1982-83. 


Figure  11.  Moisture  tension  at  Taxiway  A,  1979-80. 
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Figure  12.  Vertical  resilient  displacement  at  two  load  levels  and  three  radii.  Taxiway  B,  Test 
Point  T5. 


FROST-SUSCEPTIBILITY 
INDEX  TESTING 

Frost-susceptibility  index  tests  allow  geotechni¬ 
cal  engineers  to  determine  the  potential  for  frost 
heaving  and  thaw  weakening  of  subgrade  soils  and 
unbound  base  and  subbase  materials  in  roads  and 
airfields.  In  a  survey  of  transportation  depart¬ 
ments  throughout  the  world,  Chamberlain  et  al. 
(1984)  found  that  most  agencies  have  developed 


their  own  unique  frost-susceptibility  index  criteria 
based  on  laboratory  tests,  that  these  criteria  fail  to 
discriminate  marginally  frost-susceptible  material 
from  material  that  is  frost  susceptible,  and  that 
there  is  little  documentation  of  the  efficacy  of  the 
adopted  standards.  Furthermore,  most  of  the  vari 
ous  tests  consider  only  frost  heave  or  thaw  weak¬ 
ening  rather  than  both,  and  most  of  those  that 
employ  laboratory  freezing  tests  require  excessive 
time  and  impose  poor  control  of  test  conditions. 


n 
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The  objective  of  this  study  was  to  develop  im¬ 
proved  index  test  methods  for  fully  characterizing 
the  frost  susceptibility  of  soils.  To  accomplish  this 
task,  a  thorough  review  of  frost-susceptibility  in¬ 
dex  tests  and  practices  of  transportation  agencies 
was  made.  The  index  tests  were  categorized  into 
three  types  or  levels  of  complexity.  The  efficacy  of 
each  test  was  identified  and  the  attributes  and  de¬ 
ficiencies  were  noted. 

During  the  evaluation  phase  it  was  decided  that 
one  index  test  from  each  of  the  three  types  would 
be  selected  for  further  evaluation.  This  was  done 
so  that  a  geotechnical  engineer  could  select  an  in¬ 
dex  test  with  a  particular  level  of  reliability  and 
complexity  commensurate  with  the  size  and  scope 
of  a  project. 

The  index  tests  sheeted  include  the  U.S.  Army 
Corps  of  Engincc  s  frost  design  soil  classification 
system,  a  moisture-tension/hydraulic-conductiv¬ 
ity  test,  and  a  new  freezing  test  containing  both 
frost-heave  and  thaw-weakening  elements.  These 


tests  were  run  on  the  materials  from  the  Winchen- 
don  and  Albany  test  sites.  The  results  were  com¬ 
pared  with  field  observations  of  frost  heave  and 
thaw  weakening  at  the  test  sites,  and  the  validity 
of  each  test  was  determined. 

Details  of  the  selection  process,  the  procedures 
for  the  selected  tests,  and  the  analysis  of  the  data 
are  given  in  three  reports  by  Chamberlain  (1981a, 
1981b,  in  prep.  a). 

Index  tests  selected 

Corps  of  Engineers  method 

The  U.S.  Army  Corps  of  Engineers  frost  design 
soil  classification  system  (Berg  and  Johnson  1983) 
in  its  present  form  is  shown  in  Table  1.  It  is  based 
on  three  levels  of  screening:  (level  I)  the  percent¬ 
age  of  particles  smaller  than  0.02  mm,  (level  II) 
the  soil  type  under  the  Unified  Soil  Classification 
Sysem,  and  (level  III)  a  laboratory  freezing  test. 
The  ratings  of  frost-susceptible  soils  according  to 


Table  1.  U.S.  Army  Corps  of  Engineers  frost  design  soil  classification  system. 
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Typical  soil  type 
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PFStf 
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SI 
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this  method  are  given  in  six  categories:  negligible, 
very  low,  low,  medium,  high  and  very  high.  The 
simplest  rating  (based  on  level  1  tests)  is  the  classi¬ 
fication  of  negligible  frost  susceptibility  given  to 
gravels  with  less  than  1  .S^o  finer  than  0.02  mm 
and  sands  with  less  than  3%  finer  than  0.02  mm. 
Any  soil  failing  this  criterion  requires  comDlete 
soil  classification  tests  (level  11).  Gravels  with 
1 .5-3%  finer  than  0.02  mm  and  sands  with  3-10% 
finer  than  0.02  mm  also  require  a  laboratory  frost- 
heave  test  (level  III).  The  range  of  possible  degrees 
of  frost  susceptibility  is  very  wide  for  most  soils. 
For  this  reason,  the  Corps  of  Engineers  procedure 
includes  the  recommendation  that  a  freezing  test 
be  performed  for  more  precise  determination  of 
the  degree  of  frost  susceptibility  of  all  frost- 
susceptible  soils. 

Moisture-tension/ 
hydraulic-conductivity  test 

The  second  laboratory  test  selected  for  consid¬ 
eration  was  a  moisture-tension/hydraulic-conduc¬ 
tivity  test.  This  test  was  selected  because  it  ad¬ 
dresses  the  fundamental  causes  of  ice  segregation 
and  frost  heave  more  directly  than  particle-size 
tests  do.  It  is  also  a  test  that  is  routinely  conducted 
at  CRREL  for  research  purposes  and  was  already 
included  in  the  frost-heave  modeling  part  of  this 
research  project.  However,  the  lengthy  period  of 
time  required  to  conduct  this  test,  the  requirement 
for  a  very  skillful  technician,  and  the  inconclusive 
nature  of  the  results  led  to  the  conclusion  that  this 
test  is  not  suitable  for  determining  the  frost  sus¬ 
ceptibility  of  soils.  Details  of  this  evaluation  are 
given  by  Chamberlain  (in  prep.  a). 

New  freezing  test 

As  a  result  of  the  literature  review,  it  was  con¬ 
cluded  that  all  available  freezing  tests  have  at  least 
one  serious  fault.  The  flaws  include  too  small  or 
too  large  a  sample  size,  significant  radial  heat 
flow,  lack  of  surcharge,  insufficient  control  of 
moisture  availability,  lack  of  appropriate  frost- 
susceptibility  assessment  criteria,  and  insufficient 
field  validation.  Our  experience  with  the  C  RRIT 
freezing  test  (Kaplar  1974,  Chamberlain  and  Car- 
bee  1981)  specified  as  part  of  the  Corps  of  1  ngi 
neers  procedures  for  pavement  design  in  frost  areas 
(Berg  and  Johnson  1983)  indicates  that  it  suffers 
from  several  defects.  The  most  significant  are 
poor  temperature  control,  indeterminate  side  fric¬ 
tion,  lengthy  test  period,  lack  of  thaw  weakening 
index,  and  provision  for  only  a  single  freeze 

It  was  concluded  that  a  new  freezing  test  should 
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Figure  13.  Schematic  of  new  freezing  test  apparatus. 

be  developed,  for  which  we  established  the  follow  ¬ 
ing  basic  guidelines: 

•  The  test  should  be  as  simple  as  possible,  so 
that  highway  and  geotechnical  laboratories 
can  conduct  tests  readily. 

•  The  results  must  be  reliable. 

•  The  test  conditions  must  bear  a  relation  lo 
freezing  conditions  in  the  tield. 

•  The  test  conditions  must  also  relate  to  thaw 
weakening  in  the  field. 

•  The  test  must  be  of  short  duration 

•  The  test  must  accommodate  the  complete 
range  of  material  types  from  granular  base 
and  subbase  materials  to  fine-grained  sub 
grade  materials 

•  The  apparatus  should  be  inexpensive  to 
construct  and  operate. 

We  developed  a  new  freezing  test  using  these 
guidelines  (Fig.  13).  The  equipment  includes  a 
rubber-membrane-lined,  multi-ring  freezing  cell  to 
minimize  side  friction,  liquid-cooled  cold  plates 
for  precise  top  and  bottom  boundarv  temperature 
control,  and  a  data  acquisition  and  control  svstem 
lor  automated  temperature  control  and  data  proc 
essmg  The  tevt  imposes  two  freeze  thaw  cvcles  to 
account  lot  the  changes  in  suscepnbiluv  to  tiost 


heave  caused  by  a  prior  freeze-thaw  cycle.  Four 
samples  are  tested,  each  150  mm  in  diameter  and 
150  mm  in  height.  Samples  of  prepared  material 
are  compacted  to  field  density  and  moisture  condi¬ 
tions.  Water  is  freely  available  through  porous 
base  plates.  The  test  requires  five  days  to  com¬ 
plete.  The  heave  rate  at  the  end  of  the  first  eight 
hours  of  each  two-day  freeze-thaw  cycle  is  used  as 
an  index  of  frost-heave  susceptibility.  A  CBR  test 
is  conducted  after  the  second  thaw  to  provide  an 
index  of  thaw -weakening  susceptibility.  As  will  be 
seen,  both  indices  (frost  heave  and  thaw  weaken¬ 
ing)  must  be  used  to  determine  the  frost  suscepti¬ 
bility  of  a  soil.  Details  of  the  test  and  the  pro¬ 
cedures  are  provided  by  Chamberlain  (in  prep.  b). 

Laboratory  test  results 

Particle-size  tests,  Atterberg-limit  tests,  and  the 
new  freezing  test  were  performed  on  the  materials 
from  the  Winchendon  and  Albany  test  sites. 

Corps  of  Engineers  method 

Frost-susceptibility  ratings  according  to  the 
Corps  of  Engineers  frost  design  soil  classification 
method  are  shown  in  Table  2.  All  but  one  of  the 
soils  (Hyannis  sand)  were  rated  frost  susceptible 
by  this  method.  The  frost-susceptibility  ratings 
vary  widely  and  do  not  appear  to  be  strongly  relat¬ 
ed  to  either  the  heave  rate  or  the  pavement  deflec¬ 
tion  observed  in  the  field.  The  rating  correctly 
identifies  the  Hyannis  sand  as  non-frost-suscepti¬ 
ble  (negligible  to  low  frost  susceptibility)  in  terms 
of  the  frost  heave  measured  in  the  field,  but  the 
other  soils  require  freezing  tests  to  determine  their 
degree  of  frost  susceptibility. 


The  frost  susceptibilities  of  the  test  soils  were 
also  determined  (Table  2)  by  finding  a  similar  soil 
in  a  table  of  previous  C  RREL  freezing  test  results 
published  by  Kaplar  (1974)  and  as  part  of  Corps 
of  Engineers  guidance  on  this  subject  (Berg  and 
Johnson  1983).  The  Hyannis  sand  was  the  only 
test  soil  that  received  a  rating  of  non-frost-suscep¬ 
tible.  Others,  such  as  the  dense-graded  stone  and 
Graves  sand,  are  ranked  highly  frost  susceptible. 
Most  notable  is  the  considerably  narrower  range 
of  frost  susceptibility  for  each  soil  given  by  this 
method  than  that  obtained  by  the  particle-size  and 
soil  classification  method. 

The  frost  susceptibilities  of  the  test  soils  were 
also  determined  from  the  C'RREL  standard  freez¬ 
ing  test  in  the  tapered  cylinder.  Table  2  also  sum¬ 
marizes  the  test  results.  Again,  the  Hyannis  sand 
was  ranked  as  the  lowest  in  frost  susceptibility 
among  the  Winchendon  materials.  The  lowest  rat¬ 
ing  of  all  the  test  materials,  however,  was  for  the 
Albany  Airport  Taxi  way  B  subgrade  material 

New  freezing  test 

An  example  of  the  results  of  the  new  freezing 
test  is  shown  in  Figure  14.  Four  freezing  tests  were 
conducted  on  each  of  the  Winchendon  soils  and 
two  on  each  of  the  Albany  materials.  The  heave 
rates  (Fig.  15)  were  lowest  for  Hyannis  sand  and 
highest  for  the  Sibley  till  material.  The  high  heave 
rates  for  Sibley  till  material  occurred  during  the 
second  freeze.  Sibley  till  also  had  the  lowest  CBR 
values  after  the  second  thaw  (Fig.  16).  The  Taxi¬ 
way  B  subgrade  material  had  the  highest  post¬ 
thaw  CBR. 

Figure  15  shows  that  the  frost-heave  rales  for 


(able  2.  Summary  of  frost-susceptibility  ratings  according  to  all  criteria. 
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Figure  14.  Frost-heave  results  for  the  new  freeze-thaw  test  on  four  sum 
pies  of  Sibley  till. 
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Figure  15.  Heave  rates  for  all  test  materials  in 
ne  w  freeze -  tha  w  test. 

three  of  the  soils  were  significantly  greater  during 
the  second  freeze.  For  two  of  the  soils  (Taxiway  A 
base  and  Taxiway  B  subbase)  the  heave  rate  in¬ 
creased  by  factors  of  2  to  3,  while  for  Sibley  till  it 
increased  by  a  factor  of  9  from  one  of  the  lowest 
heave  rates  (2  mm/day)  to  the  highest  (18  mm/ 
day).  This  illustrates  the  importance  of  including 
the  second  freeze-thaw  cycle  in  this  test. 

Figure  16  shows  that  the  CBR  values  were  re- 


Figure  16.  CBR  values  for  all  test  materials  in 
new  freeze-thaw  test. 

duced  by  two  cycles  of  freezing  and  (hawing 
Again,  the  detrimental  change  was  greatest  tor 
Sibley  till,  with  the  CBR  being  reduced  to  about 
I2°?n  of  the  value  before  freezing 

Discussion 

Comparisons  ot  the  labor  itnn  and  field  (lost 
heave  rates  are  shown  m  I  mire  l~  With  the  ex 
ccption  of  the  Taxiway  \  results,  there  is  a  Miong 
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figure  17.  (  ompartson  of  laboratory 
and  field  heave  rates  during  first  and 
second  freezes. 


Figure  18.  Comparison  of  CBR  after  thaw¬ 
ing  with  maximum  resilient  pavement  deflec¬ 
tion  during  thawing. 
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correlation  between  these  heave  rates  for  the  first 
free/e  thaw  cycle.  The  correlation  is  not  on  a  line 
of  equality,  as  the  laboratory  heave  rates  exceed 
the  field  values  by  a  factor  of  10  or  more.  How¬ 
ever,  since  it  was  the  intent  of  this  study  to  use  the 
freezing  test  qualitatively,  not  as  a  quantitative 
predictor  of  frost  heave  in  the  field,  the  differ¬ 
ences  are  not  considered  significant.  When  the  re¬ 
sults  are  plotted  for  the  second  freeze,  the  correla¬ 
tion  between  the  laboratory  and  field  results  be¬ 
comes  weaker,  as  the  points  for  the  two  Albany 
base  materials  and  the  Sibley  till  fall  far  to  the 
right  of  the  curve  fitting  the  remainder  of  the  data. 

Tne  correlation  between  the  C  BR  after  thawing 
and  the  maximum  resilient  pavement  deflection 
during  thawing  (Fig.  18)  is  better  than  the  correla¬ 
tion  between  the  frost-heave  parameters.  In  this 
case,  all  the  average  values  of  deflection  fall  close 
to  a  straight  line  showing  inverse  proportionality 
with  C  BR  after  thawing. 

The  comparisons  in  Figures  \7  and  18  clearly 
show  the  need  for  including  a  thaw-weakening  in¬ 
dicator  as  a  frost-susceptibility  index  in  the  labor¬ 
atory  freezing  test.  If  only  the  heave  rate  from  the 
first  freeze  had  been  used  to  determine  the  frost 
suscep  .bility  of  the  Sibley  till,  the  soil  would  have 
been  called  non-frost -susceptible.  The  results  of 
the  second  freeze,  the  C  BR  after  thawing,  and  the 
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pavement  deflection  under  load  in  the  field  show 
that  this  soil  is  clearly  frost  susceptible. 

As  it  was  desirable  to  compare  at  least  qualita¬ 
tively  the  results  of  the  various  test  methods  for 
the  test  soils,  preliminary  frost-susceptibility  clas¬ 
sification  criteria  were  established  for  the  new 
freezing  test.  These  criteria  are  shown  in  Table  3. 
They  are  based  on  the  average  frost-heave  rate 
during  the  first  8  hours  of  the  first  freeze  and  the 
CBR  after  two  cycles  of  freezing  and  thawing.  The 
criteria  were  established  by  comparing  a  limited 
number  of  laboratory  and  field  tests,  so  they  must 
be  considered  preliminary  and  subject  to  change. 


Table  3.  Preliminary  frost-suscepti¬ 
bility  criteria  for  the  new  freezing  test. 
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Figure  19.  Comparison  of  frost-susceptibility  classifications 
by  the  methods  studied. 


r-rost-susceptibility  criteria  for  the  fieid  obser¬ 
vations  were  also  developed  for  comparison  with 
the  criteria  for  the  new  test.  Details  of  these  rat¬ 
ings  are  given  by  Chamberlain  (in  prep.  a).  Table  2 
shows  ratings  of  all  the  test  materials  using  the  cri¬ 
teria  for  the  new  freezing  tests  and  the  field  obser¬ 
vations  along  with  a  summary  of  the  ratings  using 
the  other  procedures  in  this  study.  The  bar  charts 
in  Figure  19  illustrate  the  range  of  frost-suscepti¬ 
bility  classifications  for  two  of  the  test  soils. 

The  comparisons  of  the  laboratory  results  and 
the  field  observations  show  that  the  Corps  of  En¬ 
gineers  particle-size  test  was  effective  in  rating 
non-frost-susceptible  soils  but  was  not  effective  in 
determining  the  degree  of  frost  susceptibility  of 
the  other  soils.  Estimates  made  from  published  re¬ 
sults  of  previous  CRREE  freezing  tests  predict 
frost-heaving  behavior  better.  Elowever.  it  is  clear 
from  the  test  results  for  several  of  the  study  mater¬ 
ials,  and  especially  for  the  Sibley  till  results,  that  a 


thaw-weakening  indicator  such  as  the  CBR  test 
after  freezing  and  thawing  is  required. 

The  importance  of  conducting  two  freeze-thaw 
cycles  is  not  clearly  demonstrated  in  these  tests  ex¬ 
cept  by  the  fact  that  a  prior  freeze-thaw  cycle  can 
have  a  dramatic  effect  on  the  behavior  during 
freezing  in  the  laboratory.  The  heave  rates  during 
the  second  freeze  probably  correlate  poorly  with 
the  field  heave  rates  because  prior  freeze-thaw 
cycles  at  the  field  sites  had  not  extended  through 
the  full  depth  of  the  test  material  to  the  water 
table.  The  entire  path  through  which  water  must 
How  from  the  water  table  to  the  frost  front  was 
not  preconditioned  by  prior  freezing  and  thawing 
as  it  was  in  the  laboratory.  Thus,  laboratory  frost- 
susceptibility  ratings  obtained  from  a  second 
freeze  are  probably  useful  only  when  prior  frost 
cycling  in  the  field  fully  penetrates  the  material  to 
the  water  table.  More  field  experience  is  needed  to 
fully  understand  this  problem. 
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Conclusions 

It  is  clear  from  this  study  that  to  determine  the 
frost  susceptibility  of  a  soil  accurately,  it  is  neces¬ 
sary  to  conduct  a  freezing  test.  While  the  Corps  of 
Engineers  frost  design  soil  classification  method  is 
useful  for  separating  non-frost-susceptible  soils 
from  frost-susceptible  soils,  it  is  not  very  helpful 
in  determining  their  degree  of  frost  susceptibility, 
and  it  cannot  be  effectively  used  to  predict  the  de¬ 
gree  of  thaw-weakening  susceptibility. 

The  heave  rate  in  the  laboratory  freezing  test 
can  be  used  to  determine  the  frost-heave  suscepti¬ 
bility  in  the  field,  and  the  CBR  value  after  freezing 
and  thawing  is  a  strong  indicator  of  field  thaw 
weakening  leading  to  increased  resilient  pavement 
deflection  under  load. 

The  requirement  for  two  freeze-thaw  cycles  was 
not  clearly  demonstrated  in  these  tests;  however,  it 
remains  a  candidate  element  of  the  freezing  test 
for  use  in  regions  where  freeze-thaw  cycling  ac¬ 
tively  conditions  the  soil  fabric  throughout  the  full 
depth  to  the  ground  water  table. 

The  freezing  test  proposed  is  a  feasible  candi¬ 
date  for  replacing  the  CRREL  freezing  test  as  the 
standard  for  the  Corps  of  Engineers  because  it  re¬ 
quires  much  less  time  to  conduct  (5  vs  14  days),  it 
provides  much  better  boundary  temperature  con¬ 
trol,  it  eliminates  the  side  friction  problem  preva¬ 
lent  in  the  current  standard  test,  it  provides  an  in¬ 
dicator  of  thaw-weakening  susceptibility  as  well  as 
an  indicator  of  frost-heaving  susceptibility,  and  it 
allows  the  effects  of  repeated  freezing  and  thawing 
to  be  determined. 

As  further  research  using  the  new  freezing  test 
appears  justified,  the  automated  equipment  neces¬ 
sary  to  conduct  the  test  should  be  provided  to  sev¬ 
eral  transportation  laboratories  to  obtain  a  larger 
data  base  from  which  more  reliable  frost-suscepti¬ 
bility  criteria  can  be  established.  The  Corps  of  En¬ 
gineers  frost  design  soil  classification  method 
based  on  grain  size  and  soil  classification  should 
be  included  in  any  additional  studies,  and  the  cor¬ 
responding  criteria  should  be  refined  based  on  the 
experience  obtained. 


SOIL  COLUMN  AND 
DUAL  GAMMA  SYSTEM 

Design  features 

It  has  long  been  known  that  moisture  content 
and  soil  density  change  significantly  when  soil 
freezes.  Changes  in  these  variables,  in  time  and 
space,  need  to  be  evaluated,  as  they  are  critical  pa- 
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Figure  20.  Soil  column  (schematic). 


rameters  in  heat  and  mass  transfer  during  freezing 
and  thawing.  The  soil  column  and  dual  gamma 
system  were  constructed  to  generate  such  data 
from  laboratory  tests.  The  data  were  then  used  to 
develop,  verify  and  refine  the  mathematical  model 
of  frost  heave  and  thaw  settlement. 

Three  separate  soil  columns  were  designed  and 
fabricated  for  this  study.  The  interior  dimensions 
are  about  50  mm  in  diameter,  135  mm  in  diameter, 
and  100  mm  by  100  mm  square,  respectively.  Fig¬ 
ure  20  is  a  schematic  drawing  of  the  square  soil 
column.  The  columns  consist  of  300-mm-long  see- 


Figure  21.  Dual  gamma  system. 


merits  that  are  stacked  to  provide  any  desired 
length  from  900  to  1800  mm.  The  upper  300  mm 
of  the  column  is  tapered  to  reduce  the  side  friction 
that  develops  as  freezing  occurs  from  the  top 
downward  (Berg  et  al.  1980b). 

The  base  of  each  column  is  fitted  with  a  porous 
stone  and  tubing  leading  from  a  source  of  water. 
These  features  and  a  Marriott  tube  allow  us  to 
control  the  position  of  the  water  table  in  the  col¬ 
umn.  The  soil  column  is  placed  inside  a  temper¬ 
ature-controlled  jacket  to  minimize  lateral  heat 
losses  and  promote  one-dimensional  freezing  or 
thawing. 

The  top  of  the  column  is  attached  to  an  air- 
act  is ated  piston,  which  applies  the  desired  force 
(surcharge)  to  the  soil  surface  through  the  surface 
heat-exchange  plate.  This  plate  is  connected  to  a 
refrigerated  bath  via  flexible  tubing.  The  bath 
temperature  is  set  manually  to  obtain  the  desired 
rate  of  freezing  A  thermocouple  is  embedded  in 
the  surface  heat-exchange  plate.  A  linear-motion 


potentiometer  (LMP)  is  attached  to  the  side  of  the 
soil  column  to  monitor  the  amount  of  frost  heave 
of  the  soil.  Other  instrumentation  in  the  soil  col¬ 
umn  includes  thermocouples  to  monitor  tempera¬ 
tures  at  selected  positions  and  tensiometers  to 
monitor  pore-water  pressures  (moisture  tension). 
The  amount  of  water  drawn  into  the  soil  is  moni¬ 
tored  manually  once  or  twice  per  day,  and  the  ten¬ 
siometers,  thermocouples  and  LMF’  are  attached 
to  a  data  collection  system  that  records  data  hour¬ 
ly. 

A  dual  gamma  system  was  designed  to  be  used 
with  the  soil  column.  It  nondestructive^  monitors 
changes  in  density  and  moisture  content  with  time 
during  freezing  and  thawing  of  soils  in  the  soil  col¬ 
umn.  The  system  (Cireatorex  et  al.,  in  prep.)  con¬ 
sists  of  two  nuclear  sources,  an  electronic  detector 
to  monitor  gamma  radiation  from  the  sources,  a 
tower  to  position  the  sources  and  detector  vertical¬ 
ly,  and  electronic  equipment  lo  control,  monitor 
and  record  data  from  the  system  (big  21) 
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Test  results 

Several  tests  were  run  on  soils  from  the  Win- 
chendon  site  and  other  soils,  and  the  results  were 
used  as  one  source  of  data  for  validating  and  re 
fining  the  frost-heave  model,  a  mathematical 
model  of  coupied  heat  and  moisture  flow.  During 
late  1984  and  1985,  two  special  tests  were  conduct¬ 
ed  using  the  soil  column  and  dual  gamma  system. 
Both  >pecial  tests  were  designed  to  generate  data 
to  be  used  to  validate  and  refine  the  recently  devel¬ 
oped  thaw -settlement  algorithm  of  the  frost-heave 
model. 

During  previous  tests  using  the  soil  column, 
frost  heave  was  not  substantial,  even  with  highly 
frost-susceptible  soil  in  the  column.  Therefore,  for 
these  two  tests  we  froze  150-200  mm  of  soil  in  the 
multi-ring  molds  used  in  the  laboratory  frost-sus¬ 
ceptibility  test.  Graves  sand  from  the  Winchendon 
test  sections  was  used  for  the  tests.  The  specimens 
were  150  mm  in  diameter  by  150  mm  high  and 
were  molded  at  a  density  of  about  1.52  gm/cm' 
for  the  tests.  The  soil  specimens  were  frozen  slow¬ 
ly  from  the  top  downward  in  a  temperature- 
controlled  cabinet:  during  freezing,  about  20  mm 
of  heave  occurred. 

The  frozen  soil  specimens  were  machined  to  fit 
snugly  inside  the  upper  segment  of  the  soil  col¬ 
umn,  which  was  about  135  mm  in  diameter.  The 
lower  900  mm  of  the  soil  column  was  filled  with 
unfrozen  Graves  sand.  Temperature  sensors  and 
tensiometers  were  installed  at  several  locations  in 
both  the  frozen  and  unfrozen  portions  of  the  soil 
column  An  I  MP  was  attached  to  the  surface 
heat-exchange  plate  to  monitor  settlement  as 
thawing  progressed.  A  surcharge  of  about  7  kPa 
was  applied  at  the  surface  using  the  air-actuated 
piston  The  temperature  in  the  jacket  surrounding 
the  soil  column  was  maintained  at  about  -0,5 C. 
The  water  table  was  maintained  at  a  depth  of 
about  900  mm  below  the  top  ot  the  sample 
throughout  the  tests. 

After  attaching  a  frozen  soil  segment  to  the  rest 
ot  the  soil  column  and  installing  the  instrumenta¬ 
tion.  we  circulated  a  mixture  of  ethylene  glycol 
and  water  at  3  ('  through  the  surface  heat- 
exchange  plate  for  several  davs.  This  caused  the 
frost  line  to  advance  further  and  allowed  tempera¬ 
tures  and  pore  water  pressures  in  the  soil  to 
stabilize  prior  to  the  most  important  part  of  the 
test,  whose  obiective  was  to  monitor  changes  dur¬ 
ing  thaw  for  this  purpose  the  temperature  of  the 
circulating  fluid  was  gradually  increased  to  about 
•  2  (  Thawing  progressed  xlowlv  from  the  top 
downward  during  the  tests,  taking  six  days  for 
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Figure  22.  Thaw  set  tie  men  t  in  two  soi ! -col¬ 
umn  tests  on  Graves  sand. 


complete  thawing  in  Test  1  and  about  eight  days  in 
Test  2. 

Figure  22  shows  settlement  during  thawing  in 
Tests  1  and  2,  while  Figure  2.  shows  the  changes 
in  temperature  and  moisture  stress.  Slight  positive 
pressures  wert  observed  in  both  tests  before  the 
entire  frozen  layer  was  thawed.  When  thawing  was 
complete,  water  drained  downward  into  the  un¬ 
derlying  soil.  The  results  from  the  tests  were  used 
as  data  sets  for  refining  and  validating  the  thaw- 
settlement  portion  of  the  frost-heave  model.  The 
simulated  settlements,  temperatures  and  moisture 
stresses  obtained  by  use  of  the  model  agree  very 
well  with  the  observed  values  for  Test  2.  The  simu¬ 
lation  of  Test  1  was  not  as  accurate,  mainly  due  to 
malfunctions  of  the  laboratory  surface-tempera¬ 
ture  controller  and  the  data  collection  system. 
Because  of  these  problems,  reliable  upper-boun¬ 
dary  temperature  data  were  unavailable  for  the 
simulation. 

The  dual  gamma  system  was  used  periodically 
to  monitor  the  moisture  contents  and  density  pro¬ 
files  in  the  soil  column.  Observations  were  closely 
paced  in  the  thawing  portion  of  the  column  and 
50  or  1 1 H)  mm  apart  in  the  unfrozen  portion. 


I  mure  23,  Temperatures  and  pore  pressures  m  two  soil-column  tests  on  Graves  sand.  The  solid  lines  are  meas¬ 
urements;  the  dashed  lines  are  simulations. 


MATHEMATICAL  MODEL  OF 

FROST  HEAVE  AM)  THAW  SETTLEMENT 

Model  development 

W  e  have  reported  the  development  of  the  model 
rn  Here  et  al.  (1980a,  b),  Guvrnon  et  al.  (1980, 
1981a,  1981b  and  in  prep.)  and  Hromadka  et  al. 
(1981,  1982).  I  he  model  assumes  one-dimensional 
vertical  heat  and  moisture  flux.  It  is  intended  for 
use  on  problems  of  seasonal  freezing  and  thawing 
of  nonplastic  soils  beneath  pavements  in  which 


frost  does  not  penetrate  deeply  into  soils  beneath 
the  water  table.  Furthermore,  it  is  intended  for  use 
where  surcharge  effects  are  not  large  (usually  less 
than  60  kPa). 

Assumptions 

The  main  assumptions  of  the  model  are: 

•  There  are  three  distinct  zones  in  a  freezing 
soil:  a  frozen  zone,  a  freezing  zone  and  an 
unfrozen  zone. 

•  Moisture  transport  in  the  unfrozen  zone  is 
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governed  by  the  unsaturated  flow  equation 
based  on  Darcy’s  L  aw . 

•  Moisture  flow  is  via  liquid  movement,  and 
vapor  How  is  neglible. 

•  Moisture  How  in  the  frozen  zone  is  negligi¬ 
ble,  and  there  is  no  moisture  loss  or  gain  at 
the  frozen  soil  surface. 

•  Soil  deformations  in  the  unfrozen  zone  are 
negligible. 

•  Soil  pore-water  pressures  in  the  freezing 
zone  are  governed  by  an  unfrozen-water- 
content  factor. 

•  All  processes  are  single  valued;  i.e.,  there  is 
no  hysteresis. 

•  Heat  transport  in  the  entire  soil  column  is 
governed  by  the  sensible-heat-transport 
equation  including  a  convective  term. 

•  The  frozen  zone  is  deformable  for  deter¬ 
mining  thermal  parameters;  i.e.  the  thermal 
conductivity  is  modified  to  compensate  for 
the  increased  path  length  caused  by  frost 
heave. 

•  Salt  exclusion  processes  are  negligible;  i.e.. 
the  unfrozen  water  content  is  constant  with 
respect  to  temperature. 

•  Phase-change  effects  (freezing  and  thawing) 
and  moisture  effects  can  be  modeled  as  de¬ 
coupled  processes. 

•  Freezing  and  thawing  can  be  approximated 
as  isothermal  phase-change  processes. 

•  All  parameter  and  model  uncertainty  can  be 
incorporated  into  a  universal  probability 
model  applicable  to  a  specific  class  of  soils. 

•  Fluxes  of  heat  and  water  are  vertical,  i.e. 
the  model  is  one-dimensional  and  no  lateral 
fluxes  are  considered. 

Governing  equations 

The  appropriate  equation  describing  soil  mois¬ 
ture  flow  consistent  with  the  above  assumptions  is 
derived  by  substituting  the  extended  Darcy  mois¬ 
ture-flow  law  into  the  one-dimensional  continuity 
t  luation  for  an  incompressible  fluid  in  porous 
media: 

ft  (ttt  ()j  rWj 

—  ! A|>  Hh/Hx)  =  +  -  —  (I) 

where  the  total  hydraulic  head  h  equals  the  sum  of 
the  pore  pressure  head  (hp  -  u  i>w,  where  u  is  the 
pore-water  pressure  and  is  the  density  of  w  ater) 


side  of  the  equation,  o,  is  the  density  of  ice,  is 
the  volumetric  ice  content,  and  «u  is  the  volumetric 
water  content.  The  ice  sink  term,  iyW,/ OvO'T  ex¬ 
ists  only  for  a  freezing  or  thawing  zone,  and  in 
these  zones  equation  1  is  coupled  to  the  heat-trans¬ 
port  equation.  The  ice-sink  term  assumes  that  ",  is 
a  continuous  function  of  time. 

It  is  computationally  convenient  to  represent 
the  soil-water  characteristics  as  a  known  or  as¬ 
sumed  function  relating  pore-water  pressure  u  and 
volumetric  water  content  tiu.  This  is  done  by  deter¬ 
mining  point  values  of  0U  and  u  in  the  laboratory 
and  by  least-squares  fitting  of  a  function  similar  to 
that  used  by  Gardner  (1958)  to  the  data: 


u  -4*  hp  "  +  1 


where  n  is  the  initial  porosity  and  a  and  A*  are 
best-fit  parameters  for  each  soil. 

Similarly  the  eoefficient-of-permeability  func¬ 
tion  for  an  unsaturated  soil  is  determined  in  the 
laboratory  by  determining  point  values  of  and 
u  for  that  soil  and  by  least-squares  fitting  of  Gard¬ 
ner’s  function  to  the  data: 


A'm  -  ~ 


and  the  elevation  head  (/ic 


v).  The  coefficient 


of  permeability  A,,  is  a  function  ol  the  pore  pres¬ 
sure  head  in  the  unfrozen  soil  zone.  In  the  right 


Hk  hp  '  +  1 


where  A's  is  the  saturated  permeability  and  AK  and 
d  are  the  best-fit  parameters. 

A  phenomenological  relationship  is  assumed  for 
adjusting  the  coefficient  of  permeability  for  un¬ 
frozen  soils  to  represent  conditions  where  ice  is 
partially  blocking  soil  pores  (in  the  freezing  zone). 
It  is  assumed  that  the  coefficient  of  permeability 
in  the  freezing  and  frozen  zones  is 

A'f  =  A'H(Ap)  •  10  E"'  when  £«,  >  0  (4) 

where  £  is  a  parameter  determined  empirically  . 
Taylor  and  L,uthin  (1978)  and  Jame  (1978)  used 
similar  concepts  to  reduce  the  permeability  in  the 
freezing  zone. 

Frost  heave  is  estimated  fror.,  the  total  amount 
of  ice  segregation  in  the  frozen  zone  by 

<K  -  «i-(n-«r)  (5) 

where  -  volumetric  segregated  ice  content 
n  -  porosity 

--  residual  unfrozen  water  content. 

If  'A  ■  0,  ice  segregation  has  occurred  and  the 


T  =  Temperature 

Tl  =  Lower-boundary  temperature 
Ty  =  Upper  boundary  temperature 
h  =  Total  hydraulic  head 
hP  =  Pore-water  tension 
.i0  =  Overburden  pressure 
«u  =  Water  content  in  the  unfrozen  zone 
»,  =  Ice  content 

",  -  Residual  unfrozen  water  content 
»5  =  Segregated  ice  content 
hPL  =  Pore  water  tension  at  the  lower  boundary 
n  =  Porosity 
D  =  Element  length 
x  =  Vertical  coordinate 
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Figure  24.  Typical  model  simulation  result  at  a  given  time. 


frost  heave  is  computed  by  multiplying  by  the 
zone  thickness.  Thaw  consolidation  from  ice  melt¬ 
ing  is  the  reverse  process. 

Figure  24  illustrates  the  solution  of  a  freezing 
problem  during  a  time  step.  The  0r  parameter  es¬ 
tablishes  the  pore-water  pressure  at  the  freezing 
front  for  the  solution  of  the  moisture  transport 
equation.  The  lower  hydraulic  boundary  condi¬ 
tion  is  usually  the  water  table.  Surcharge  and  over¬ 
burden  effects  ct0,  which  tend  to  restrain  frost 
heave,  are  modeled  by  reducing  the  pore-water 
tension  at  the  ice-segregation  front  by  the  pressure 
of  the  overlying  material  and  surcharge. 

Conceptual  basis  for 
thaw-settlement  algorithm 

The  thaw -settlement  portion  of  the  model  is  dis¬ 
cussed  separately  because  of  the  importance  of 
this  submodel  for  the  evaluation  of  thaw  weaken¬ 
ing  of  soils  beneath  pavements,  a  major  objective 
of  this  project.  The  concepts  advanced  by  Mor- 
gcnstern  and  Nixon  ( 1971 )  provide  the  framework 
for  the  thaw-settlement  and  pore-water-pressurc 
algorithm  presented  here.  In  the  past  the  limited 
availability  and  poor  quality  of  laboratory  data  in¬ 
hibited  the  development  of  accurate  and  tested 
thaw  -settlement  models.  Additional  data  were  col¬ 
lected  during  this  study  using  the  CRRFI  soil  col¬ 


umn  to  test  the  adopted  thaw-settlement  algo¬ 
rithm. 

The  Morgenstern  and  Nixon  model  is  based  on 
well-known  theories  of  heat  conduction  and  linear 
consolidation  of  compressible  soils.  Terzaghi’s 
one-dimensional  consolidation  theory  was  applied 
to  develop  a  moving-boundary  solution  applicable 
to  permafrost  soils  that  thaw  and  consolidate 
under  the  application  of  load.  A  closed-form  solu¬ 
tion  was  obtained. 

Our  application  of  the  model  is  directed  to  lay¬ 
ered  systems  of  well-compacted  soils.  Consequent¬ 
ly  the  application  is  restricted  to  winter  heaving  of 
subgrade  soils  and  spring  thaw  settlement  origin¬ 
ating  in  those  same  soils  with  no  net  consolidation 
or  change  in  pavement  elevation  occurring  over  a 
sequence  of  several  years  of  freeze-thaw  action. 

A  departure  from  the  Morgenstern  and  Nixon 
model  is  the  solution  of  the  linear  governing  equa¬ 
tion  of  excess  pore-water  pressure  (Terzaghi’s 
equation)  numerically  rather  than  exactly.  The  nu¬ 
merical  code  for  this  solution  already  exists  in  the 
frost-heave  model.  This  method  allows  more  flexi¬ 
bility  in  handling  the  upper-surface  pore-water- 
pressure  boundary  condition.  A  second  departure 
from  the  Morgenstern  and  Nixon  method  is  the 
use  of  a  more  general  heat  transport  equation. 

lo  determine  the  thaw  settlement,  equation  I. 


which  governs  the  moisture  flow  during  frost 
heave,  was  modified  to  include  a  temporal  void- 
ratio  term  (Lambe  and  Whitman  1979): 


■'ll  „  dh  t?i  dv\ 

IT  +Sm^ai  +  ^  ~aT 


where  S  is  the  degree  of  saturation  and  is  the 
coefficient  of  volume  compressibility. 

Equation  6  is  the  basis  of  the  algorithm  for  esti¬ 
mating  thaw  settlement  and  pore-water  pressure 
during  and  after  thawing.  When  soil-surface  tem¬ 
peratures  are  above  freezing  and  the  upper  ele¬ 
ment  is  fully  saturated,  soil-surface  pore-water 
pressures  are  set  to  a  specified  pressure,  which  is 
usually  atmospheric  plus  the  excess  pressure  from 
overburden  effects.  However,  the  model  does 
have  the  ability  to  apply  a  specified  lesser  positive 
pressure  representing  drainage  upward  through  a 
slowly  leaking  pavement  overlying  the  base  course 
material.  When  the  upper  soil  element  becomes 
partly  drained  (5  <  100c/o)  or  when  the  surface 
element  refreezes,  the  soil-surface  boundary  con¬ 
dition  for  the  moisture  equation  is  reset  to  a  no¬ 
flux  boundary  condition. 

As  thawing  progresses  downward,  each  discrete 
soil  element  is  checked  to  determine  the  degree  of 
saturation.  If  excess  pore-water  pressure  exists, 
w  ater  in  excess  of  the  initial  porosity  is  treated  as  a 
source,  forcing  an  upward  drainage  of  water.  The 
hydraulic  conductivities  of  underlying  frozen 
zones  are  determined  by  equation  3,  with  the  pa¬ 
rameters  Ak  and  J  evaluated  by  laboratory  tests 
and  the  pore-water  pressure  determined  by  means 
of  equation  2  with  0U  set  equal  to  6r. 

When  the  soil  column  is  completely  thawed  and 
reconsolidated,  free  downward  drainage  occurs  in 
accordance  with  equation  1,  and  the  no-flux  sur¬ 
face  boundary  condition  is  assumed  to  persist. 

Boundary  and  initial  conditions 

The  model  requires  auxiliary  conditions  as  fol¬ 
lows: 

•  Initial  conditions  for  pore  pressure,  ice  con¬ 
tent  and  temperature. 

•  Soil-surface  boundary  conditions  for  pore 
pressure  and  temperatures  (which  may  vary 
with  time). 

•  l  ower  boundary  conditions  for  pressure  and 
temperature  (which  may  vary  with  time). 

While  there  are  many  possibilities  tor  incorpotal 
ing  boundary  conditions  in  the  model  to  suit  spe 
cific  applications,  the  CRREL  version  of  the 
model  has  the  features  discussed  below. 


The  upper-boundary  pore-water  pressure  is 
fixed  at  a  specific  value  with  respect  to  time  during 
freezing.  Prior  to  freezing  it  varies  according  to 
equation  1,  while  after  thawing  it  varies  according 
to  equation  6. 

The  lower  pore-pressure  boundary  condition  is 
usually  a  set  of  discrete  pore-water  pressures  (ten¬ 
sions)  at  specific  times  that  are  related  to  the  water 
table  elevation.  Intermediate  times  and  pore-water 
pressures  are  linearly  interpolated. 

The  upper  temperature  boundary  condition  con¬ 
sists  of  a  set  of  specified  step  functions  such  as 
mean  daily  air  temperatures.  These  values  are 
multiplied  by  a  factor  to  represent  soil-surface 
temperatures,  in  a  manner  similar  to  the  Corps  of 
Engineers  n-factor  approach  for  seasonal  freezing 
indexes. 

Bottom  temperature  boundary  conditions  con¬ 
sist  of  a  set  of  times  and  temperatures.  Tempera¬ 
tures  are  linearly  interpolated  at  intermediate 


Numerical  approach 

Numerical  solution  of  the  governing  equations, 
subject  to  their  respective  boundary  and  initial 
conditions,  is  by  the  nodal  domain  integration 
method  (Hrornadka  et  al.  1982).  The  one-dimen¬ 
sional  solution  domain  is  divided  into  a  number  of 
variable-length  finite  elements  where  parameters 
are  assumed  to  be  temporarily  constant  but  may 
vary  from  element  to  element  and  from  time  to 
time.  Figure  25  illustrates  the  division  of  a  vertical 
column  into  elements  and  nodes.  The  state  varia- 
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l  iguie  25.  Nonuntjonn  sod  profile  divided  into 
elements  and  nodes  and  showing  boundary  enndt 


ble  in  each  element  is  assumed  to  be  described  by  a 
linear  basis  function  such  that  the  state  variable  is 
continuous  throughout  the  solution  domain.  The 
time-domain  solution  is  by  either  the  Crank-Nicoi- 
son  approach  or  the  fully  implicit  method. 

Computation  is  initiated  by  giving  the  initial 
conditions,  and  the  solution  is  advanced  in  time. 
At  specified  times  nonlinear  parameters  are  up¬ 
dated. 

Probabilistic  concepts 

If  this  one-dimensional  model  were  perfect,  it 
would  predict  exactly  the  amount  of  frost  heave 
and  frost  penetration  at  a  point  with  time.  How¬ 
ever,  it  would  not  account  for  spatial  variations  in 
the  amounts  of  frost  heave  and  the  degree  of  thaw 
weakening  caused  by  differences  in  the  thermal 
and  hydraulic  properties,  in  turn  caused  by  differ¬ 
ences  in  density,  moisture  content  and  soil-particle 
size  and  arrangement.  We  assume  that  the  model 
is  not  perfect  and  combine  uncertainties  in  the 
model  with  variations  in  the  soil  via  probabilistic 
methods.  The  uncertainty  is  arbitrarily  grouped 
into  four  general  areas: 

•  Errors  due  to  the  choice  of  the  model,  in¬ 
cluding  the  choice  of  a  numerical  analog. 

•  Errors  due  to  spatial  and  temporal  varia¬ 
tions  in  the  soil  properties. 

•  Errors  due  to  inadequate  boundary  condi¬ 
tions  and  to  the  choice  of  initial  conditions. 

•  Errors  due  to  the  selection  of  parameter  val¬ 
ues. 

Errors  due  to  the  choice  of  a  model  are  proba¬ 
bly  indeterminate  by  strictly  analytical  methods. 
Choosing  a  model  is  probably  best  left  to  experi¬ 
ence  with  application  of  the  model  for  a  number 
of  known  laboratory  and  field  conditions.  Errors 
associated  with  the  choice  of  a  numerical  analog 
were  readily  examined,  as  were  errors  associated 
with  spatial  and  temporal  variations  in  the  soil 
properties  (Guymon  et  al.  1981b  and  in  prep., 
Hromadka  et  al.  1982).  We  concluded  that  one 
numerical  technique  had  little  advantage  over 
others  and  the  model  is  relatively  insensitive  to 
spatial  discretization.  The  model  is  very  sensitive 
to  temporal  changes,  and  we  found  that  for  most 
solutions  a  time  step  of  0.2  hr  and  a  parameter  up¬ 
date  frequency  of  1  hr  are  satisfactory. 

Errors  associated  with  inadequate  boundary 
and  initial  conditions,  and  particularly  with  pa¬ 
rameter  values,  require  special  attention  due  to  the 
probabilistic  nature  of  their  variations.  For  this 
reason,  probabilistic  theory  is  applied  to  consider 
the  uncertainties  and  variations. 


Most  investigations  of  this  nature  use  hundreds 
or  even  thousands  of  computer  simulations.  This 
type  of  stochastic  analysis  can  be  very  expensive, 
particularly  if  variations  arc  significantly  different 
for  different  soil  types  or  boundary  or  initial  con¬ 
ditions.  The  probabilistic  method  used  in  the 
model,  discussed  in  detail  in  Guymon  et  ai.  (1981b 
and  in  prep.),  eliminates  the  need  for  such  large 
numbers  of  computer  simulations. 

Model  verification 

Model  verification  has  been  a  continuing  proc¬ 
ess  since  completion  of  early  work  on  formulating 
the  model  reported  by  Berg  et  al.  (1980a).  This  re¬ 
port  contained  early  verification  of  decoupled 
components  of  the  model  (e.g.  sensible  heat  trans¬ 
port)  against  analytical  solutions  using  linearized 
computer  simulations  (Fig.  26).  As  verification 
work  progressed,  it  was  found  necessary  to  refine 
the  computer  code  to  more  accurately  simulate 
pore  pressures,  temperatures  and  frost  heave 
Guymon  et  al.  (1980)  further  reported  on  verifica¬ 
tion  efforts  using  Fairbanks  silt  as  a  test  case.  Sub¬ 
sequently  Guymon  et  al.  (1981a,  1981b  and  in 
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Figure  26.  Measured  and  simulated  frost 
depth  and  frost  heave,  laboratory  soil-column 
lest  on  Fairbanks  silt. 
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Figure  2  7.  Measured  and  simulated  frost  depth  and  frost  heave.  Taxiway  B,  A  Ibany  County 
Airport,  1979-80. 


prep.)  presented  much  greater  detail  in  verifying  tent  are  selected  so  that  the  moisture  tension  is 

the  model  against  laboratory  and  field  data.  Fig-  75-100  kPa.  The  actual  values  may  be  soil  specific 

ure  27  illustrates  calculated  and  measured  values  and  much  greater  or  lower  than  this  range, 

of  frost  heave  and  frost  penetration  with  time  Even  if  more-precise  scientific  knowledge  were 

from  one  of  the  test  sections  at  the  Albany  County  available  for  this  function,  calibration  might  still 

Airport.  Additional  verification  results  are  pre-  be  required.  There  is  no  model  in  existence  for 

sented  below  .  porous-media  flow  processes  that  does  not  require 

calibration  to  improve  the  confidence  limits  so 
Discussion  that  the  deterministic  solution  becomes  relatively 


The  model  produced  good  results  for  soils  rang¬ 
ing  from  silts  to  relatively  coarse-grained  and  mar¬ 
ginally  frost-susceptible  soils.  Moreover,  these  re¬ 
sults  have  been  compared  with  carefully  controlled 
laboratory  data  as  well  as  less-precise  field  data 
for  three  locations. 

To  achieve  such  results,  however,  good  esti¬ 
mates  of  hydraulic  parameters  are  required.  Judg¬ 
ment  is  required  in  assigning  appropriate  values, 
as  there  may  be  considerable  variations  in  even  the 
most  carefully  measured  soil  parameters,  particu¬ 
larly  unsaturated  hydraulic  conductivity. 

The  modeling  exercise  also  requires  calibration 
of  the  hydraulic-conductivity  correction  factor  E, 
a  phenomenological  correction  factor  for  soil  in 
the  freezing  zone.  This  factor  is  used  as  the  pri¬ 
mary  calibration  parameter  to  achieve  the  results 
presented  here.  If  measurements  of  hydraulic  con¬ 
ductivity  for  partly  frozen  soil  were  available  or 
could  be  easily  obtained,  it  is  likely  that  those  data 
could  be  used  instead  of  the  phenomenological  re¬ 
lationship  now  incorporated  in  the  model. 

The  model  also  requires  an  estimate  of  the  mois¬ 
ture  tension  in  the  freezing  zone.  This  is  done  indi¬ 
rectly  by  selecting  a  residual  (unfrozen)  water  con¬ 
tent  for  the  frozen  zone  and  by  calculating  the  cor¬ 
responding  pore  water  tension  by  the  Gardner 
equation.  Generally  values  of  residual  water  con- 


precise.  Hypothetical  solutions  of  such  problems 
using  assumed  parameters  have  a  considerable  er¬ 
ror,  which  for  some  porous-media  problems  may 
be  tolerable  in  the  engineering  analysis  process. 
Usually  engineering  judgment  and  experience  are 
exercised  to  infer  the  level  of  certainty  of  such 
computations.  This  need  is  evident  in  the  problem 
considered  here. 

We  believe  that  the  model  simulates  phenomena 
in  the  freezing  zone  adequately  for  our  present  en¬ 
gineering  purposes  and  that  it  will  meet  the  need 
of  practicing  pavement  engineers  for  predicting 
frost  heave  and  some  of  the  parameters  influenc¬ 
ing  thaw  weakening  of  pavement  systems.  The  de¬ 
velopment  of  a  model  more  closely  linked  to  ac¬ 
cepted  concepts  of  soil  physics  awaits  a  more  com¬ 
plete  understanding  and  formulation  of  processes 
in  the  freezing  zone.  It  would  also  require  addi¬ 
tional  computer  time  and  expense  to  solve  a  more 
complex  formulation  of  the  processes,  and  addi¬ 
tional  time,  equipment  and  expense  for  conduct¬ 
ing  laboratory  tests  to  define  additional  soil  pa¬ 
rameters. 

The  output  from  the  model  includes  cumulative 
frost  heave  with  time  at  the  surface,  subsurface 
temperatures  and  pore-water  pressures.  The  pre¬ 
dicted  frost  heave  can  be  used  directly  to  aid  in  se¬ 
lecting  an  appropriate  pavement  design  by  relating 


it  to  pavement  roughness  criteria.  Temperatures 
are  used  to  determine  the  positions  of  the  freezing 
and  thawing  zones,  and  temperatures  and  pore- 
water  pressures  are  used  in  empirical  equations  de¬ 
veloped  from  laboratory  tests  to  estimate  resilient- 
modulus  values  for  layers  within  the  pavement 
system  at  various  times  of  th.  .ar.  The  resilient- 
modulus  data  are  then  used  in  a  pavement  struc¬ 
tural  response  model,  where  output  can  be  related 
to  pavement  performance  criteria.  Examples  of 
this  application  are  presented  later  in  this  report. 

SEASONAL  VARIATION  IN  THE 
RESILIENT  MODULUS  OF 
GRANULAR  SOILS 

In  areas  of  seasonal  frost  the  supporting  capaci¬ 
ty  of  subgrade  soils  and  unbound  base  and  sub¬ 
base  materials  for  roads  and  airfields  can  vary 
widely  during  freeze-thaw  cycling  and  the  subse¬ 
quent  spring  and  summer  recovery  period.  In  this 
phase  of  the  research  we  sought  to  evaluate  the 
seasonal  fluctuations  in  material  properties  for  six 
soils  from  the  Winchendon  test  site  and  five  soils 
from  pavements  at  the  Albany  County  Airport 
test  site.  The  investigations  focused  on  the  under¬ 
lying  cause  of  premature  distress  in  pavement  sys¬ 
tems  that  are  susceptible  to  frost— the  reduction 
of  the  resilient  modulus  of  subgrade  soils  and  un¬ 
bound  base  courses  during  and  following  spring 
thaws.  In  this  context  the  resilient  modulus  is  con¬ 
ventionally  defined  as  the  deviator  stress  divided 
by  the  resilient  strain  (i.e.  the  recoverable  strain). 
The  research  was  concerned  with  frost-susceptible 
granular  soils  exhibiting  little  or  no  cohesion  and  a 
high  degree  of  nonlinear  (i.e.  stress-dependent) 
mechanical  behavior.  The  research  objective  was 
to  develop  laboratory  methods  for  characterizing 
the  seasonal  changes  in  the  resilient  modulus  of 
such  materials  throughout  a  complete  annual  cy¬ 
cle.  Field  tests  were  conducted  to  validate  the  lab¬ 
oratory  methods. 

Repeated-load  triaxial  tests  were  performed  to 
determine  the  resilient  characteristics  of  the  mater¬ 
ials  under  conditions  simulating  those  during  the 
field  tests.  The  laboratory  triaxial  tests  were  per¬ 
formed  on  soils  in  the  frozen,  thawed,  recovering 
and  recovered  conditions.  Empirical  relationships 
were  then  generated  by  standard  statistical  tech¬ 
niques  to  express  the  resilient  modulus  Mr  as  a 
function  of  the  density,  the  soil-moisture  tension 
and  the  stresses  imposed  in  the  triaxial  tests.  For 
frozen  soil  and  asphalt  concrete,  temperature  is 
also  a  key  parameter. 


For  validation  purposes,  field  tests  were  used  to 
determine  the  surface  deflection  response  ol  paved 
soil  test  sections  under  plate  loads.  Surface  deflec¬ 
tion  basins  were  measured  under  loads  imposed  h\ 
a  repeated-load  plate-bearing  (RPB)  apparatus 
and  a  falling-weight  deflectometer  (FWD)  1  he 
tests  were  performed  at  critical  times  between  late 
fall  and  late  spring  to  characterize  the  variation  in 
load  response  throughout  the  freeze-thaw -recov¬ 
ery  cycle. 

The  validity  of  the  laboratory  results  was  then 
examined  by  comparing  the  measured  deflection 
basins  with  deflection  basins  calculated  tor  the  test 
section  using  the  expression,  for  resilient  modulus 
developed  from  the  laboratory  tests.  In  using  these 
expressions,  temperatures  and  moisture  tensions 
measured  at  the  time  of  each  field  loading  test 
were  applied  to  evaluate  the  resilient  modulus 
Layered  elastic  analyses  of  the  test  sections  under 
the  conditions  prevailing  during  each  field  loading 
test  yielded  stresses,  strains  and  resilient  vertical 
displacements  throughout  the  system;  calculated 
surface  deflection  basins  were  thus  generated  and 
compared  to  the  deflection  basins  actually  meas¬ 
ured  in  the  field. 

Interim  results  were  given  by  Cole  et  al.  (1981) 
and  Johnson  et  al.  (1982).  Detailed  procedures,  re¬ 
sults  and  analyses  of  repeated-load  triaxial  tests, 
and  of  field  in-situ  plate-loading  tests  and  the  cor¬ 
responding  deflection  basin  analyses,  are  given  in 
a  four-part  report  series  by  Cole  et  al.  (1986  and  in 
prep.)  and  Johnson  et  al.  (1986a,  b). 

Characterization  by  laboratory  testing 

The  objective  of  this  phase  of  the  work  was  to 
develop  procedures  for  obtaining  realistic  expres¬ 
sions  for  the  resilient  modulus  of  granular  soils  in 
the  frozen,  thawed  and  recovered  states  in  terms 
of  the  significant  variables.  The  testing  and  ana 
lytical  techniques  we  have  developed  allow  us  to 
simulate  in  the  laboratory  the  gradual  recovery  of 
stiffness  experienced  in  the  field  as  a  frost-suscep¬ 
tible  soil  drains,  consolidates  and  desaturates  after 
thawing. 

Experimental  approach 

The  experimental  approach  called  for  repeated 
load  triaxial  tests  on  all  the  asphalt  concrete  and 
test  soils  from  both  test  sites.  We  obtained  field 
cores  of  the  asphalt  concrete  and  of  the  finer- 
grained  soils  in  the  frozen  state.  We  also  sampled 
several  tests  sections  at  the  Winchendon  site  be 
fore  the  subsoils  froze,  to  characterize  the  material 
after  it  had  fully  recovered  from  the  pievious  win 
ter’s  freezing  cycle.  In  the  case  of  glacial  till  or 


gravel,  where  core  samples  could  not  be  obtained, 
we  reconstituted  specimens  in  the  laboratory  and 
subjected  them  to  open-system  freezing  at  a  rate 
of  25  mm  day  to  produce  the  appropriate  condi¬ 
tions  for  testing 

1  oad  cycles  on  the  soils  were  applied  using  two 
waveforms  to  simulate  the  loading  pulses  associ¬ 
ated  with  our  two  field-testing  devices:  the  falling- 
weight  deflectometer  (FWD)  and  the  repeated- 
load  plate-bearing  (RPB)  apparatus.  Two  hundred 
load  cycles  were  applied  at  each  level  of  confining 
and  deviator  stress  used  in  the  laboratory  tests, 
and  the  resilient  modulus  and  resilient  Poisson’s 
ratio  were  calculated  when  a  nominally  steady- 
state  response  was  achieved. 

Bv  performing  these  tests  under  a  variety  of 
temperature  and  moisture  conditions  comparable 
to  those  observed  in  the  field,  we  were  able  to 
characterize  the  resilient  properties  and  generate 
analytical  expressions  in  terms  of  the  significant 
variables. 

The  terminology  used  to  describe  the  state  of 
the  soils  is  as  follows.  “Frozen”  refers  to  material 
with  at  least  some  pore  ice  present,  and  "thawed” 
refers  to  material  having  only  liquid  pore  water 
but  that  is  still  suffering  from  the  effects  of  a 
freezing  cycle  b  e.  it  has  not  drained  or  consolidat¬ 
ed  to  its  prefree/ing  state  and  is  consequently  still 
in  a  weakened  condition).  The  term  “recovered” 
refers  to  material  t hai  by  means  of  drainage,  con¬ 
solidation  and  desaturation  has  been  restored  to 
the  moisture  and  density  onditions  (and  hence 
stillness)  prevailing  before  the  start  of  a  Ireezing 
cycle  Details  of  the  experimental  procedures  are 
given  in  Cole  el  al.  ( 1986) 

Advancements 

A  number  of  significant  advancements  in  triax- 
lal  testing  equipment  and  procedures  have  allowed 
us  to  simulate  the  thaw  and  recovery  process  in  the 
laboratory  I  he  kev  variable  in  the  recovery  proc¬ 
ess  is  soil -moisture  tension.  As  a  soil  drains  after 
thawing,  it  first  reconsolidates  to  a  condition  of 
zero  pore  water  pressure  Following  this  phase, 
gradual  desaturation  occurs,  the  moisture  tension 
rises,  and  the  resilient  modulus  increases  We  de 
veloped  a  svstem,  using  removable  inaxial -cell 
bases  equipped  with  tens.ometers,  that  allowed  us 
to  retest  a  given  specimen  several  times  at  increas¬ 
ing  levels  of  moisture  tension,  thus  simulating  the 
changes  observed  in  the  field  I  ach  specimen  re 
mamed  mounted  on  its  base  throughout  the  test 
me  sequence,  and  excessive  handling  was  thus 
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higure  28.  Resilient  modulus  vs  unfro¬ 
zen  water  content  for  three  frozen 
H’inchendon  soils  tested  at  69  kPa  con¬ 
fining  stress  and  the  indicated  ranges 
o  f  deviator  stresses. 

The  ability  to  perform  sequential  tests  allowed 
us  to  simulate  the  stiffness  recovery  in  the  labora 
tory  with  a  relatively  small  number  of  specimens. 
We  generally  tested  specimens  four  times,  which 
allowed  us  to  cover  the  full  range  of  field  condi¬ 
tions  with  each  specimen.  The  use  of  moisture  ten¬ 
sion  as  the  primary  means  of  describing  the  soil 
state  has  proven  to  be  effective  since  it  strongly  in¬ 
fluences  the  resilient  modulus  and  is  relatively  eas¬ 
ily  monitored  in  both  the  field  and  the  laboratory. 

Other  developments  at  CRRFI  in  the  use  of  the 
pulsed  nuclear  magnetic  resonance  technique  on 
frozen  soils*  provided  us  with  information  on  the 
unfrozen  water  content  of  the  frozen  soils  tested  in 
this  program  The  unfrozen  water  content  (Fig 
28)  was  found  to  have  a  profound  influence  on  the 
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Figure  29.  Resilient  modulus  vs  tern 
perature  for  the  si.x  Winchendon  test 
sods  m  the  frozen  state.  Curves  are  gener 
a  led  hy  regression  equations  based  on  un¬ 
frozen-  water-content  expressions. 

resilient  modulus  of  material  in  the  frozen  state 
and  was  consequently  used  as  the  key  variable  in 
the  empirical  modeling  effort.  The  equations  are 
of  the  form 

V/r  =  C'  (wu/w,)c  (7) 

where  Vfr  =  resilient  modulus 
C  ,C:  =  constants 

wu  -  unfrozen  water  content,  expressed 
a>  a  function  of  temperature 
w,  =  total  gravimetric  water  content. 

The  unfrozen  water  content  is  expressed  (Cole 
1984)  hy  an  equation  of  the  form 

wu  -  .»<  T)  (8) 

where  o  and  are  regression  constants  and  T  is 
temperature  This  constitutes  the  first  use  of  an 
untrozen-water  content  term  in  a  resilient-modu¬ 
lus  expression,  and  the  results  have  been  very 


satisfactory.  The  expressions  are  relatively  simple 
and  mathematically  well  behaved,  and  they  can  be 
extrapolated  to  lower  temperatures  with  reasona 
hie  confidence.  This  analytical  approach  requires 
fewer  laboratory  tests  to  obtain  a  useful  expres¬ 
sion.  Figure  29  shows  a  plot  of  modulus  vs  tem¬ 
perature  for  the  six  Winchendon  soils. 

Standard  statistical  techniques  have  been  used 
in  this  work  to  generate  empirical  expressions  for 
the  resilient  modulus  in  terms  of  the  unfro/en 
water  content  for  the  frozen  case  and  in  terms  of 
the  moisture  tension,  the  applied  stress  and  in 
some  cases  the  dry  density  for  the  thawed  cases. 

In  some  cases  we  have  applied  the  commonly 
used  bulk-stress  model  for  the  stress  dependency 
of  a  nonlinear  material,  which  is  of  the  form 

W,  =  K.Jk  (9) 

where  A,  and  A.  are  constants  and  J  is  the  first 
stress  invariant  (sum  of  the  principal  stresses, 
equal  to  the  bulk  stress  ft).  We  have  also  employed 
a  somewhat  more  complex  stress  function  involv¬ 
ing  the  second  stress  invariant  and  the  oetahcsi  .u 
shear  stress,  of  the  form 

M,  =  A,(/./rocl)A  (10) 

where  J.  =  second  stress  invariant 

—  0,n,  +  (!.<>,  f  '1(1, 

Toct  =  octahedral  shear  stress 

=  ''«(((), +  (<>:-<l,);  4  (”  <■,)-] 

=  principal  stresses. 

The  stress  function  of  equation  10  is  unique  !ia 
it  accounts  for  the  effects  of  both  confining  pres¬ 
sure  and  the  principal  stress  ratio  on  the  modulus 
in  a  manner  appropriate  for  many  granular  mater¬ 
ials. 

Figure  30  shows  a  plot  of  modulus  vs  bulk  stress 
for  a  thawed  Winchendon  test  soil  (Hyannix  sand). 
The  bulk-stress  model  does  not  account  lot  the 
fact  that,  for  certain  types  of  soil,  the  modulus  de¬ 
creases  with  increasing  principal  stress  ratio.  The 
stress  function  given  in  equation  10.  however,  ac¬ 
counts  for  the  influence  of  the  stress  ratio  and  thus 
gives  a  more  efficient  linear  representation  ol  the 
data  (Fig.  31). 

Results 

Analysis  of  the  laboratory  test  data  tor  each  soil 
yielded  equations  applicable  to  the  frozen  state 
and  the  thawed  state.  The  frozen-state  equation  is 
in  fact  valid  up  to  and  including  the  point  i*  com- 
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Figure  30.  Resilient  modulus  of  thawed  Hyannis  sand  vs  J, 
(first  stress  invariant,  equal  to  bulk  stress  0)  for  several  values 
of  principal  stress  ratio,  demonstrating  the  tendency  for  the 
modulus  to  decrease  with  increasing  stress  ratio  when  the  con¬ 
fining  stress  o,  remains  fixed. 
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f  igure  31.  Resilient  modulus  vs  J:  /roc,  for  the  test  data  given  in 
f  igure  30. 


plete  thaw ,  at  which  point  the  expressions  yield  the 
average  modulus  value  observed  in  the  “as 
thawed”  condition  (i.e.  prior  to  any  drainage).  A 
significant  feature  of  this  approach  is  that  it  yields 
a  continuous  f.nction  up  to  the  point  of  thawing, 
because  we  included  appropriate  data  points  from 
the  thawed  state  in  the  analysis. 

The  laboratory  testing  sequence  mentioned 
above  generated  data  that  allowed  us  to  model  the 
post-thaw  recovery  process  with  a  single  equation 
employing  terms  for  stress,  soil-moisture  tension 
and  dry  density.  The  increasing  stiffness  associ¬ 
ated  with  the  recovery  phase  is  predicted  through 


the  moisture-tension  term,  which  is  incorporated 
in  the  coefficient  K,.  Figure  32  shows  how  mois¬ 
ture  tension  affects  the  term  K,,  and  thus  the  mod¬ 
uli,  in  equations  9  and  10  as  evaluated  by  means  of 
regression  analyses  (Tables  4  and  5)  for  several  Al¬ 
bany  Airport  test  soils.  Earlier  work  (Cole  et  al. 
1981)  demonstrated  that  K,  is  statistically  indepen¬ 
dent  of  the  moisture  tension. 

Tables  4  and  5  show  the  results  of  the  regression 
analyses  on  the  asphalt  concrete  and  al)  test  soils 
for  both  the  frozen  and  thawed  conditions.  In  gen¬ 
eral,  the  stress  function  J>/t0 c,  produced  a  higher 
coefficient  of  determinatioi.  (R1)  than  the  J,  (bulk 
stress)  model. 


Table  4.  Results  of  regression  analyses,  asphalt  concrete  and  test  soils  from  Winchendon. 


Material 

Load  pulse 

Regression  equation 

Asphalt  concrete 

RPB 

Af.(MPa)  =  exp[9. 204-5. 552  x  10  T-9.744  *  10  *P] 

Haversine 

A/,(MPa)  =  exp(9. 183-7.4? >.  lO-’rif®1’'' 

FWD 

A/r(MPa)  =  exp(9. 429-7. 47x  lOT] 

Natural  subgrade 

RPB  and  FWD 

AT,(MPa)  =  8.829  f,(af  '<* 

RPB  and  FWD 

ATr(MPa)  =  20.74  /;<»)°  152 

Grates  sand 

Frozen 

RPB 

VTr(MPa)  =  exp(9. 677-1 .0314T-0 .0708  T'M  rK/  o0)'° 

RPB 

A/, (MPa)  =  39.1(w/w,)  1 

FWD 

V/r<MPa)  =  32.14(wu/w,)  1  “ 

Thawed 

RPB 

A/,(MPa)  =  2.l39x  10‘/(v) ; '',2’ /.(o)U4W 

FWD 

A/r(MPa)  =  9.27  x  10’  /(A)  ;  60  /, ( <z)°  477 

RPB 

Wr(MPa)  =  6.68  x  10*  /(v)  J  •**' /;(u)°414 

FW'D 

A/,(MPa)  =  1 .47  x  10*  /(v  >  1 ''  /:(o)04" 

Recov  ered 

RPB 

Vfr(MPa)  =  6.89 /(<j)041* 

RPB 

AT,(MPa)  =  4.80  /;(,j)‘,4U4<’ 

Ikalanian  sand 

Frozen 

RPB 

Af,(GPa)  =  exp(13.74-(0.820)r-(0.0538)r-(0.8378)w 

i-(0.04416)w;](r<s,/,Jo)  0  )l! 

RPB 

A/r(MPa)  =  86.4(tvu/w,) ■' 

Thawed 

RPB 

A/,(MPa)  =  8.l29x  10‘/(v-)  !  ,:4/(-,d>"  r,/i(<r)04W 

RPB 

Mr(MPa)  =  3.02!  xl04/(.;)  >!“/(yd)"6,4/,(or)0"2 

Recovered 

RPB 

AT, (MPa)  =  5.69  x  104/(A)'!  IIS  /,(o)°  557 

RPB 

M, (MPa)  =  2.405  x  I04/(v)  j*"/.(o)°“! 

Hart  Brothers  sand 

Frozen 

FWD 

AT,(MPa)  =  38.28(wa/w,)  1 ’,J 

RPB 

ATf(MPa)  =  4.085  x  10'  (w„/w,>  1  ” 

FWD 

.AT, (MPa)  =  8.05x  10- /bd)’ *7 

FWD 

Af,(MPa)  =  4.689 x  10-  (wjw,)  '  >• 

Thawed 

RPB 

AT,(MPa)  =  2.97  x  10'  /(A)'J  /(7d)?  ”6  f:(of 4S) 

RPB 

A/r(MPa)  =  1 .269  x  10’  /{A)'!  0,5  /(yd)’ 07 1  fi(of 451 

FWD 

Mr(MPa)  =  3  .93  x  104/(AV2  67  /(yd)6  "  f\of 

FWD 

M,{ MPa)  =  3.81  x  104  MY2  “7  /(yd)’ 4) 

Hyannis  sand 

Frozen 

RPB 

Af,(MPa)  =  0,68  /(yd)"  0  (wu/w,)  !  12 

RPB 

ATr(MPa)  =  33  45(wu/w[)-'0> 

Thawed 

FWD 

ATr(MPa)  =  7, 147  x  I04  fur' 7,2  f,(°f  M 

FWD 

A/,(MPa)  =  3.J7x  10’/('0',r‘/.(<J)<”*“ 

Dense-graded  stone 

Frozen 

RPB 

■Afr(MPa)  =  82.27(wu/w,)i01 

Thawed 

RPB 

A/,(MPa)  =  1.56x10  '  MY' '*  f,(of 

RPB 

A/,(MPa)  =  7,I7x  10*  yXv) •'  fAof 

Sibley  till 

Frozen 

RPB 

ATr(MPa)  =  1  01  x  lO'fw/w,) 

Thawed 

RPB 

Af,(MPa)  =  7  47,  10*  /(v)2,;',/.-('00l,; 

RPB 

M,(MPa)  =  1  29  x  10'  f{i)  1  M 

NOTES: 

RPB  repeated-load  plate-bearing  apparatus  waveform 

FWD  =  falling-weight  deflectometer  waveform 

n  =  number  of  points 

/:(,T)  ’(nil  *ol 

R ■  a  coefficient  of  determination  V  =  l 

V/f  =  resilient  modulus 

-  unfrozen  water  content 

f  =  load  \ 

*ave  frequency 

w.  ^  total  water  content 

/(  >  -  U 

,) 

T  -  44 

•  -  stress 

rkPa) 

••  -  temperature  (  (.  \ 

85 


Sul 

t.q 

R: 

error 

no 

0.97 

0.287 

1 

0  81 

0  469 

V 

3 

0.67 

0  235 

4 

0.76 

0.201 

s 

0.88 

0  332 

7 

0  91 

0  502 

8 

0.95 

0.446 

9 

0.76 

0.209 

10 

0.71 

0.224 

1  1 

0  89 

0.144 

12 

0.86 

0.15? 

13 

0.-6 

0.24? 

14 

0.87 

0  185 

15 

0,90 

0.308 

16 

0.92 

0.'49 

1  7 

0.84 

(I  '23 

18 

0.89 

0.2 '6 

19 

0  88 

0.205 

20 

0.84 

0.238 

21 

0  95 

0.53 

22 

0.92 

0.623 

n 

0.97 

0.445 

23 

0  96 

0.464 

25 

0.71 

0.280 

26 

0.87 

0.185 

27 

0.67 

0.292 

28 

0.67 

0.292 

29 

0.96 

0.536 

30 

0.95 

0.61" 

11 

0  ?| 

0.129 

12 

0.74 

0  194 

33 

0.97 

0.4)3 

34 

0  65 

0.202 

15 

0  65 

0.203 

36 

0  8? 

0.71 

.17 

0.63 

0.283 

38 

0  54 

0.313 

19 

Kim  '8  .1  .  ,| 
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Table  5.  Results  of  regression  analyses,  asphalt  concrete  and  test  soils  from  Albany  Airport. 


S(d. 

Eq 

Material  Load  pulse 

Regression  equation 

n 

R '■ 

error 

no. 

Taxiwa?  A 

Asphalt  concrete  FWD 

+Afr(  M  Pa) 

= 

1.84s  10-  exp(-3  80 x  10  ;r-9.l4x  10  T) 

88 

0.97 

0.19 

1 

RPB 

A/(MPa) 

= 

1 .01  a  10*  exp(-6.50 <  10  ;r-6.50x  10‘T] 

93 

0.98 

0.24 

s 

Haversine 

Afr(MPa) 

= 

1,09. 10*  exp(-4.75x  10  :r-7.81  x  10  T:l/S;° 

280 

0.97 

0.22 

3 

Thawed  base  FW  D  RPB 

+A/,(MPa) 

= 

UOx  10*l/(v>|  -  “  /:(t’)° 

222 

0.82 

0.16 

4 

Af.'MPa) 

= 

4.44  a 

222 

0.82 

0.16 

5 

A/, (MPa) 

3.68  .x  10*[/(v)I  :,s/:(*>)0H’/('m)344 

222 

0.84 

0.16 

6 

A/r(MPa) 

= 

2.56  «  IO*[/(v)|  ,'”/1(.-)0,'/(>d);® 

222 

0,82 

0.16 

7 

Frozen  base 

+A/,(MPa) 

1.89x  lO'lw^/w,)  4*;,  w’u  =  3  x  1  O' -( -  T)  °  ’5,  wt  =  0.075 

78 

0.78 

0.66 

8 

Thawed  Aubbase  FWD  RPB 

+A/,(MPa) 

= 

2.07x101/(0]  ’‘"/(o)029 

149 

0.80 

0.20 

9 

A/r(MPa) 

= 

4.35  x  10“l/(c)|  2  2  /(<>)° JT 

149 

0.80 

0.20 

10 

Af,(MPa) 

= 

1.39  x  10  “!/(■:  )|- ■’ u  /,(■>)"  29/<-,d) ' 00 

149 

0,82 

0  20 

1 1 

Afr(MPa) 

- 

h.oo  a  io*L/(C» '■  /(-!dr'  ” 

149 

0.82 

0.19 

12 

Frozen  subbase 

+A/r(  MPa) 

= 

8.18  x  10  (wu.  w,)  402,  wu  =  3x  10  ;(-n  o:s.  w,  =  0.055  53 

0.70 

0.84 

13 

Non-trozen  FWD  RPB 

fA/,(MPa) 

1  34  x  I0*l/(c))  1  ‘V^)0  '5 

262 

0.80 

0.80 

14 

■■ubgrade 

Taxiwax  B 

A/r(MPa) 

= 

'.73xlO’l/(c)) ■' 54 /•(<>)“” 

262 

0.78 

0.17 

15 

Thawed  base  FW  D  RPB 

*A/r(MPa) 

= 

5.55  x  10  »(/(v)] 4 

173 

0.69 

0.26 

16 

subba.se 

A/r(MPa) 

= 

9.67  x  10*(/(v)|  4  ,\/(u)ow' 

173 

0.73 

0.24 

17 

A/r(MPa) 

= 

4.28 x  I0‘l/(C)]  r/(ii),y ’ 

173 

0.71 

0.25 

18 

A/(MPa) 

= 

1 .56  x  104(/(O|  ,w/(o)°>‘/(7d)'  '2 

173 

0.74 

0.23 

19 

Frozen  base  subbase 

+  Vf,(  MPa) 

= 

l.OOx  lO'IV1*!)  iM.  WU  =  3x  10  :(-Do::,  w,  =  0.05 

92 

0.96 

0.42 

20 

Thawed  subgrade  FWD  RPB 

+A/r(MPa) 

= 

8.76  x  10’  L/(v)l  2  '*  MM0  )0 

293 

0.72 

0.20 

21 

A/(MPa) 

= 

3.36x10*  L/(01  :  15 /,(a)°  54 

293 

0.68 

0.21 

22 

Vfr(MPa) 

= 

3.80  x  104(/(v)|  * 16 /i(e)'5  25  /(>dl  * 06 

293 

0.74 

0.19 

23 

A/r(MPa) 

= 

1 .35  x  10*[/(d)]'2 13  /(o)°  ,4/(7,j)'j  06 

293 

0.70 

0.20 

24 

Frozen  subgrade 

A/(MPa) 

= 

2.66(wu/w wu  =  3.14xlO:(-7y  0  29. 
w>,  =  0.29 

152 

0.82 

0.92 

25 

V/(MPa) 

= 

2 . 59( wu / w,)'° ®5/, ( o)°  91 .  =  3.14x  10:(-7y0  29, 

w,  =  0.29 

152 

0.84 

0.85 

26 

A/,(MPa) 

= 

3.31  x  10'(wu/w,)  0  87/,(ci)06*,  =  3.14x  10-'(-7T°29. 

2.66<wu/w,)-,0V;(o-)0’8.  wu  =  3.14x  10 '(-7T029. 

152 

0.82 

0.92 

27 

Frozen  subgrade 

AT,  (MPa) 

= 

152 

0.82 

0.92 

25 

w,  =  0.29 

,W,(MPa) 

= 

2.59<wu/wt)-°*!/.(o)09\  wu  =  3. 14  x  10  '(-7”)  0  29. 
w,  =  0.29 

152 

0.84 

0.85 

26 

A/r(MPa) 

= 

3.31  x  lO'fWu/w,)  0*7/,(o)068,  =  3.14x  10 ‘(-r/0  29. 

w,  =  0.29 

152 

0.82 

0.92 

27 

Nonfrozen  subgrade 

A/(MPa) 

= 

5.16x  10*  [/(Ol  2'l/I(O026 

278 

0.81 

0.15 

28 

A/(MPa) 

= 

5.48  x  I0‘[/(C)I'2  71  /.(o)026 

278 

0.72 

0.18 

29 

Wr(MPa) 

- 

2.49  x  10*l/(v)|-2 25  U°? 26  /(',d)20’ 

278 

0.82 

0.14 

30 

N.OTPS 

RPB  -  repeated-load  pb.e-oeanng  apparatus  waveform 

f  W D  =  falling-weight  detlectomeier  waveform 

f  =  equations  i.sed  tn  analysis 

/(.)  -  <101. 36-0 J  = 

first  stress 

invariant  (kPa) 

n  =  number  of  points 

.  -  moisiure  lension  (kPa)  J,  = 

second  stress  invariant  ( k Pa ) 

R  coefficient  of  determination 

...  »  1  kPa  'o.,  = 

octahedral  shear  stress  (kPa) 

Vf.  '  thenl  modulus 

fin)  =  (J  /(-, j)  = 

■'  d  '■  o 

J 

/■(■')  =  (/ -d  = 

dry  unit  weight  (Mg 

m‘) 

••  temperature  (  (  ) 

/(  •)  -  = 

l  Mg  m' 

,4  |  (  .1  =  stress  (kPa)  h'u  =  unfrozen  water  content  (decimal) 

j.  load  waveform  frequents  (Hz)  -  1  kPa  H’i  =  total  water  content  (decimal) 
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Figure  32.  Coefficient  K,  of  the  empirical  models  (eq  9  and  Hi)  v.v 
moisture  tension  for  several  soils  from  Taxi  way  B,  Albany  Airport. 


Field  verification 

Field  in-situ  data  for  the  test  sections  from 
which  the  samples  had  been  obtained  were  needed 
for  two  purposes.  On  the  one  hand  it  was  desired 
to  measure  surface  deflections  under  plate  loads  to 
compare  with  deflections  calculated  with  the  help 
of  the  laboratory-derived  expressions  for  resilient 
modulus.  This  comparison  could  verify  the  validi¬ 
ty  of  the  laboratory  testing  methods.  Field  data 
were  also  needed  to  document  the  seasonal  varia¬ 
tion  of  temperature  and  moisture  tension  through¬ 
out  a  complete  annual  cycle.  These  variables  sig¬ 
nificantly  affect  the  resilient  modulus,  and  their 
definition  is  essential  for  a  time-dependent  evalua¬ 
tion  of  the  resilient  modulus,  which  laboratory 
testing  alone  cannot  provide. 

Deflection  measurements 

For  measuring  load-induced  deflections,  two 
types  of  in-situ  tests  were  employed  (Johnson  et 
al.  1982).  The  first  was  a  repeated-load  plate-bear¬ 
ing  (RPB)  test.  The  equipment  is  mounted  in  the 
center  of  a  large,  enclosed  semi-trailer.  In  our  tests 
the  ioad  actuator  applied  pressures  in  the  range 
from  about  200  to  600  kPa  through  a  304-mm-di- 
ameter  plate.  The  pulse  duration  was  about  1  s 
and  the  cycle  time  was  about  3  s.  Pavement  deflec¬ 
tion  was  monitored  at  radial  distances  by  means  of 
LVDTs  mounted  on  a  reference  beam. 

The  second  type  of  in-situ  test  equipment  was 
the  falling-weight  deflectometer  (FWD).  With  this 
device  a  28 -ms  pulsed  load  was  applied  to  a  300-m- 
diameter  plate  resting  on  the  pavement  surface. 
This  load  puUe  simulated  a  truck  wheel  movng  at 
moderate  peed.  The  sticks  imposed  ranged  from 
200  to  I  VitO  kPa.  I  he  pavement  deflection  at  ra¬ 


dial  distances  from  the  load  plate  was  monitored 
by  means  of  velocity  transducers. 

The  equipment  described  above  was  used  at  the 
Winchendon  test  sections  on  22  occasions  from 
October  1978  through  April  1980.  At  Albany 
County  Airport  the  same  equipment  was  used  on 
12  occasions  between  November  1978  and  June 
1980.  Between  November  1982  and  May  1983, 
FWD  tests  were  performed  11  times  at  Albany 
County  Airport. 

The  vertical  resilient  displacements  obtained 
from  the  tests  were  plotted  as  deflection  basins 
(Fig.  33).  The  displacements  at  the  Winchendon 
test  sections  measured  first  in  the  autumn  of  1978 
(Fig.  5)  decreased  to  small  amounts  in  the  second 
series  of  tests,  made  in  February  1979  when  the 
test  soils  were  frozen.  The  plots  show  the  sharp  in¬ 
crease  in  displacement  after  thawing  started.  The 
increase  in  surface  deflection  upon  thawing  was 
particularly  great  for  the  four  test  soils  at  Win¬ 
chendon  containing  the  greatest  fractions  of  fines 
(material  passing  the  No.  200  sieve),  the  Ikalan- 
ian.  Graves  and  Hyannis  sands  and  the  Sibley  till. 
Even  in  these  soils  a  substantial  decrease  in  deflec¬ 
tion  (recovery)  was  observed  within  10-20  days 
after  thawing  started.  In  the  Hart  Brothers  sand 
and  the  dense-graded  stone  the  increased  deflec¬ 
tion  upon  thawing  did  not  differ  greatly  from  the 
autumn  (recovered)  deflection.  The  comparative 
response  of  the  six  soils  to  thawing  is  illustrated  in 
Figure  34. 

Similar  plots  for  the  pavements  at  Albans 
County  Airport  also  show  the  sharp  increase  in 
displacement  after  thawing  starts  dig.  12)  At 
Taxiwav  A.  however,  the  displacement  Jut  me 
thaw  nig  was  not  much  higher  than  that  icsordcd 
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Future  33.  Deflection  basins  measured  in  1980  FH'D  test. 
Hurt  Brothers  sand  lest  section. 
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the  previous  November,  and  it  tended  to  inerease 
on  the  dates  succeeding  the  thaw ,  probably  reflect¬ 
ing  a  stronger  dependence  on  the  temperature  of 
the  thick  asphalt  concrete  pavement  than  on  the 
condition  of  the  base  or  subbase. 

Analytical  approach 

The  objective  in  analyzing  the  plate-bearing 
tests  was  to  calculate  the  deflection  basin  resulting 
from  the  application  of  a  known  load,  using  for 
each  layer  the  linear  or  nonlinear  resilient  modulus 
and  Poisson’s  ratio  determined  in  advance  by  lab¬ 
oratory  tests.  Comparison  of  the  calculated  de¬ 
flection  basins  with  the  basins  actually  observed  in 
field  tests  permitted  laboratory  test  procedures  to 
be  validated.  The  computer  program  for  analyzing 
nonlinear  elastic  layered  systems  (NELAPAV) 
chosen  for  the  calculation  of  deflections  was  de¬ 
veloped  by  Irwin  (Irwin  and  Johnson  1981)  as  a 
modification  of  the  CHEVRON  code. 

NELAPAV  incorporates  five  basic  models  of 
resilient-modulus  stress  dependency.  Any  one  of 
these  models  could  be  selected  for  any  pavement 
layer.  The  models  include 

Type  1:  V/r  =  constant 

Type  2:  A/r  =  K  JK 

Type  3:  Afr  =  A':  -HA'.-ud)A',,  nd  <  K 

A/r  =  A:  +  tod-A  )A,,  (7d  >  A i 

Type  4:  Mx  =  K  U,  roct)A 

Type  5:  M,  =  A',(roc()* 

where  A  ,  A'.-,  A,  and  A,  are  regression  constants. 
In  analyzing  the  six  Winchendon  test  sections,  we 
used  the  type  1  model  for  the  asphalt  concrete, 
types  1, 4  and  5  for  the  frozen  test  soils,  and  types 
2  and  4  for  the  thawed  and  recovered  test  soils  and 
the  natural  subgrade  test  soil.  As  the  analyses  pro¬ 
ceeded,  the  type  4  model  became  preferred  over 
the  type  2  model  for  thawed  soils.  In  analyzing  the 
deflection  on  Taxiways  A  and  B,  the  type  I  model 
was  used  for  asphalt  concrete  and  for  frozen  base, 
while  the  type  J  model  was  used  in  all  other  cases. 

In  addition  to  selecting  a  model  for  the  resilient 
modulus  and  inputting  the  applicable  values  for 
the  regression  constants,  we  selected  appropriate 
values  for  the  resilient  Poisson's  ratio  from  experi¬ 
ence  with  other  materials  and  from  published  test 
data  IT  able  6). 


Another  important  parameter  that  the  -et 
must  select  for  each  layer  is  the  coefficient  of  lat¬ 
eral  earth  pressure  at  rest  A,,  which  is  used  in  cal¬ 
culating  the  stresses  generated  by  the  weights  of 
the  lasers  overlying  the  point  for  which  the  calcu¬ 
lations  are  being  made.  Values  chosen  for  A  at 
the  Winchendon  test  sections  were  1.5  for  the  as¬ 
phalt  concrete  and  subgrade  and  1.0  for  the  test 
soils.  Tor  the  analysis  of  deflections  at  Taxiways  A 
and  B,  to  preclude  instability  of  the  siiess  func¬ 
tion,  we  chi  -e  A  values  of  1 .0  for  the  asphalt  con¬ 
crete  and  subgrade  and  0.”  for  the  base  and  sub¬ 
base 

Calculated  deflections  compared 
with  measurements 

Lor  the  deflection  calculations  the  temperatures 
and  moisture  tensions  prevailing  in  the  pavements 
at  the  time  of  each  test,  together  with  the  plate 
pressures  measured  for  each  test,  were  taken  as 
given  values,  and  the  stress-dependency  model  ap¬ 
propriate  to  each  layer  was  selected.  NELAPAV 
calculated  stresses,  strains,  displacements  and 
stress-compatible  moduli  throughout  the  system 

The  comparisons  between  the  calculated  and 
measured  deflections  at  Winchendon  for  the  high¬ 
est  plate  loading  are  summarized  in  Figure  35  for 
three  of  the  test  sections;  the  comparison  is  similar 
for  the  other  test  sections  and  for  the  lower  loads 
(Johnson  et  al.  1986a).  Several  general  trends  are 
apparent  in  these  plots.  The  maximum  deflections 
at  the  basin’s  center,  calculated  by  NELAPAV, 
tend  to  agree  well  with  the  maximum  surface  de¬ 
flections  measured  in  both  the  RPB  and  FWD 
tests,  but  they  did  not  agree  as  well  at  the  outer 
radii.  The  agreement  is  somewhat  better  for  FWD 
tests  (1989  data)  than  for  RPB  tests  (1978-1979 
data).  Perhaps  the  most  significant  observation  is 

Table  6.  V  alues  of  Poisson’s  ratio  used  in  analysis. 
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figure  35.  Measured  surface  deflections  compared  with  deflections  calculated  by  NPLAPA  I  for  three  test 


sections. 
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f  igure  36.  Measured  surface  deflections  (hWD  data)  compared  with  deflections  calculated  fn  V/  /  .1/ 
1982-83. 
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/  idlin'  3  7.  Interpretation  of  seasonal  variation  in  the  resilient 
modulus  of  six  test  soils  directly  beneath  asphalt  pavement  under 
200-  to  3 00- k  Pa  plate  pressure. 


the  reasonably  good  agreement  ot  the  post-thaw 
basin>  in  general.  The  calculated  and  measured 
deflections  differ  more  when  the  cross  section  in¬ 
cluded  lasers  of  frozen  soil.  This  problem  can  be 
attributed  in  part  to  uncertainties  in  defining  the 
exact  thickness  of  the  frozen  layers. 

The  comparisons  between  the  calculated  and 
measured  deflections  (FWD  tests)  for  Taxiway  A 
test  point  N2  and  Taxiway  B  test  point  T5  are 
shown  irt  Figure  36.  Corresponding  plots  for  the 
other  test  points  of  Taxiway  A  are  given  by  John¬ 
son  el  al  ( 1 986b).  Agreement  is  excellent. 

The  calculated  resilient  moduli  and  other  results 
trom  the  analysis  of  the  test  sections  are  summar¬ 
ized  by  Johnson  et  al.  (1986a,  b).  The  resilient 
moduli  of  the  test  soils  calculated  by  NELAPAV 
show  the  expected  seasonal  variation:  extremely 
high  values  in  the  frozen  condition,  decreasing 
dramatically  upon  thawing  and  increasing  some¬ 
what  during  the  late  spring,  summer  and  fall.  An 
interpretation  of  the  variation  in  modulus  ot  the 
upper  layer  of  each  test  soil  under  the  lower  of  the 
test  loads  is  given  in  F  igure  37  for  the  Winchendon 
test  sections.  Figures  38  and  39  show  similar  inter¬ 
pretations  tor  t tie  various  layers  in  Taxiways  A 
and  B  at  Albany  County  Airport. 

I  he  agreement  ot  the  calculated  deflections 
with  the  deflections  measured  under  F  WD  loading 
is  strong  evidence  th;  e  equations  for  nonlinear 
resilient  modulus  dov.  ed  trom  laboratory  tri- 
axial  tests  represent  v  did  characterizations  of  the 
materials  m  the  layered  pavement  system.  Accord¬ 


ingly  the  procedures  for  the  laboratory  repeated- 
load  triaxial  tests,  including  testing  at  successive 
levels  of  moisture  tension  to  track  the  recovery 
process  in  thaw-weakened  soils,  are  considered  to 
be  valid.  Those  procedures  should  be  useful  for 
structural  evaluation  and  design  of  pavements  af¬ 
fected  by  freezing  and  thawing. 

This  investigation  has  demonstrated  that  the  re¬ 
silient  modulus  depends  strongly  on  the  tempera¬ 
ture  and  moisture  tension.  These  seasonally  vary¬ 
ing  parameters  can  be  evaluated  by  installing  sen¬ 
sors  at  various  depths  below  a  paved  surface  and 
collecting  data  over  a  complete  annual  cycle.  Al¬ 
ternatively  the  frost-heave  model  of  Guymon  et  al. 
(in  prep.)  can  be  used  to  predict  both  temperatures 
and  moisture  tension  as  variables  in  time  and 


Summary  of  predictive  approav  h 

This  research  has  defined  and  verified  an  ap¬ 
proach  for  evaluating  the  resilient  modulus  of 
granular  subgrade  and  base  course  materials  as  a 
variable  in  time  and  space.  While  the  principal  fo¬ 
cus  of  the  work  was  seasonal  variation,  the  space 
variable,  represented  by  the  position  of  a  material 
in  the  pavement  profile,  also  proved  to  be  signifi¬ 
cant,  not  only  because  seasonal  variation  in  tem¬ 
perature  and  moisture  tension  itself  depends  on 
depth,  but  because  depth  below  the  pavement  sur¬ 
face  also  is  a  determinant  of  stress,  which  in  a 
granular  material  strongly  governs  its  resilient 
modulus.  Pavement  design  and  evaluation  are  af- 
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triaxial  tests  at  various  levels  of  confining  Calculations  were  made  with  the  one-dimen- 

pressure  and  deviator  stress  equivalent  to  sional  mathematical  model  to  simulate  the  frost 

the  range  of  stresses  expected  within  the  heave  and  thaw  consolidation  at  the  field  test  sites, 

pavement.  Monitor  axial  and  radial  recov-  The  time-dependent  temperatures  and  pore  pres- 

erable  and  permanent  strain,  as  well  as  sures  calculated  for  nodal  points  within  the  pave- 

stress  levels.  ment  substructure  were  used  to  determine  unique 

•Successively  desaturate  the  specimens  layers  within  the  freezing  or  thawing  system  and  to 

through  a  range  of  levels  of  moisture  ten-  select  equations  characterizing  the  resilient  modu- 

sion,  conducting  repeated-load  tests  at  each  lus  of  each  layer.  The  appropriate  equations  for 

level.  Measure  the  moisture  tension  before  resilient  modulus  were  then  input  into  the 

each  test.  NELAPAV  program  along  with  temperatures  and 

•  Through  multiple  linear  regression  tech-  pore  pressures  (moisture  tension),  and  the  defor- 

niques,  develop  equations  for  characteriz-  mations  and  stresses  at  each  point  of  interest  were 

ing  the  resilient  modulus  in  terms  of  an  ap-  calculated.  The  test  of  the  efficacy  of  the  proce- 

propriate  stress  function;  moisture  tension,  dure  required  a  comparison  of  the  observed  and 

and  possible  dry  density,  should  be  included  calculated  frost  heave,  and  the  seasonal  variation 

as  independent  variables.  of  pavement  deflection  under  repeated  loading. 

•  Evaluate  the  seasonal  variation  of  moisture 

tension  in  the  existing  or  planned  pavement  Method  of  evaluation 

substructure.  The  evaluation  can  be  made  The  procedure  for  calculating  the  frost  heave 

by  implanting  tensiometers  beneath  the  ex-  and  pavement  deflection  was  as  follows.  For  a 

isting  pavement,  or  a  similar  pavement  near-  particular  test  section,  the  appropriate  boundary 

by,  and  monitoring  them  throughout  one  conditions  and  material  properties  were  input  into 

full  year.  Alternatively  the  frost  heave  the  frost-heave  model.  Repeated  simulations  with 

model  can  be  used  to  predict  the  variation  varying  values  of  the  correction  factor  £  were  con- 

in  moisture  tension  throughout  the  thawing  ducted  with  the  mathematical  model  until  the  cal- 

and  recovery  period.  culated  frost  heave  values  agreed  with  the  field 

•  The  characterization  of  the  seasonal  varia-  values.  The  frost  depth  and  frost  heave  (or  thaw 

tion  in  the  modulus  is  essentially  complete  depth  and  thaw  settlement)  and  the  temperature 

with  the  development  of  the  regression  and  pore-water  pressure  (or  tension)  were  then 

equations  outlined  above  and  the  evalua-  printed  out  on  a  daily  basis.  The  soil  profile  for 

tion  of  the  moisture  tension  as  a  time-  each  day  was  divided  into  frozen  and  thawed 


dependent  variable.  The  resilient  modulus 
at  a  particular  time  and  under  particular 
traffic  loading  conditions  can  be  evaluated 
numerically  as  part  of  the  analysis  of  the 
pavement  cross  section.  A  suitable  pave¬ 
ment  response  model  would  be  used  to  cal¬ 
culate  the  stresses  and  the  stress-compatible 
moduli,  as  well  as  other  outputs  such  as 
strain  and  displacement. 


SEMI  EATING  FROST  HEAVE 
AM)  PAVEMENT  DEFLECTION 


The  mathematical  model  for  frost  heave  and 
thaw  consolidation  and  the  nonlinear  layered  elas¬ 
tic  pavement  response  model  NELAPAV,  in  con¬ 
cert  with  the  results  of  laboratory  resilient-modu¬ 
lus  tests,  were  used  to  simulate  frost  heave  and 
pavement  deflection  in  the  field.  Comparison  with 
field  observations  provided  the  ultimate  test  of  the 
modeling  and  laboratory  testing  procedures. 


zones.  The  thawed  zone  was  further  divided  into 
layers  according  to  moisture  tension,  and  the  fro¬ 
zen  zone  was  divided  into  layers  according  to  tem¬ 
perature.  Because  the  moisture  tension  or  temper¬ 
ature  varied  within  each  layer,  average  values  were 
calculated.  Next,  the  appropriate  expressions  for 
the  temperature-dependent  or  moisture-stress- 
dependent  resilient  modulus  A/r  were  selected  for 
the  frozen  and  thawed  soil  and  the  asphalt  con¬ 
crete  from  the  array  of  equations  given  in  Tables  4 
and  5. 

The  appropriate  models  of  stress  dependency 
contained  within  NELAPAV  were  selected,  the  re¬ 
gression  coefficients  and  other  layer  data  were  in¬ 
put  into  the  NELAPAV  program  for  each  layer, 
and  the  pavement  deflection  under  a  260-kPa  ap¬ 
plied  stress  was  calculated  for  each  day  or  time  in¬ 
terval  during  which  significant  changes  occurred 
within  the  pavement  test  section. 


Results  and  discussion 

Figure  40  shows  the  results  for  the  Graves  sand 
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Figure  40.  Comparison  of  simulated  frost  heave,  frost  and 
thaw  depths,  and  pavement  deflection  with  field  observations. 
Graves  sand  test  section,  1978-79. 


test  section  at  the  Winchendon  site.  The  calculated 
pavement  deflections  agree  well  with  the  measured 
values.  The  predictions  made  to  date  with  this  pro¬ 
cedure  demonstrate  the  efficacy  of  the  frost-heave 
and  thaw-consolidation  model  and  the  nonlinear 
layered  elastic  pav  ement  response  model  for  calcu¬ 
lations  of  seasonal  frost  heave  and  pavement  de¬ 
flection  under  load.  The  models  thus  provide  an 
instrument  for  implementing  much  of  the  work 
performed  under  this  multi-phase  research  pro¬ 
ject. 


It  is  still  necessary,  however,  to  adjust  the  E  fac¬ 
tor  in  the  mathematical  model  of  frost  heave  to 
make  the  calculated  and  observed  values  of  frost 
heave  agree.  We  expect  that  after  more  experience 
is  gained  with  the  frost-heave  model  a  correlation 
between  E  and  a  soil  factor  such  as  particle  st/e 
will  be  found  and  that  eventually  E  values  will  be 
available  as  a  chart  or  nomograph. 

Once  the  appropriate  E  factor  is  selected,  there 
appear  to  be  no  more  significant  limitations  on  the 
use  of  this  method  to  predict  frost  heave,  thaw  set 
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t lenient  and  pus  email  deflection  under  wheel 
loading  \  di-arnct  advantage  o!  this  approach  is 
that  the  tmal  results  make  it  possible  to  applv 
modern  pavement  design  sv stems  including  mech¬ 
anistic  analvsrs  and  a  cumulative  damage  ap 
pi  oach 

SI  MMAR\  OF  HNDIM.S 

froxl-suxceptibiliiy  index  tests 

\  comprehensive  review  of  the  existing  frost 
siisceptibilitv  index  tests  led  to  the  selection  ot  the 
three  most  primtising  methods  Aider  Further  eval¬ 
uation.  development  ot  new  methods,  testing  and 
comparison  with  Field  observations,  two  methods 
of  differing  levels  of  amiplexity  were  found  to  be 
useful  as  indicators  of  Field  perlormance  The 
C  orps  ot  1  ngineers  Frost  design  soil  classification 
method,  comprising  very  simple  fest  procedures  is 
useful  tor  separating  non-Frost-susceptible  from 
Frost  susceptible  soils,  but  it  does  not  reliable 
identify  their  degrees  ot  susceptibility  to  Frost 
heave  nor  o  it  effective  tor  predicting  the  degree 
ot  thaw  weakening  susceptibility.  At  a  much  high¬ 
er  level  ot  complexity,  a  newlv  developed  Freezing 
test  can  be  used  to  determine  the  Frost-heave  sus- 
ceptibilitv  in  the  Field,  and  the  C  BR  value  after 
thawing  is  a  strong  indicator  of  Field  thaw  weaken¬ 
ing  leading  to  increased  resilient  deflection  under 
applied  loads 

>uil  column  and  dual  gamma  syslem 

A  soil  column  and  dual  gamma  system  were  de¬ 
veloped  tor  use  as  research  instruments  For  obtain¬ 
ing  data  on  the  changes  in  soil  moisture  content 
and  density  that  occur  during  Freezing  and  thaw¬ 
ing  I  he  results  From  tests  on  soils  From  Winchen- 
don  and  other  sites  were  used  for  improving  and 
validating  the  Frost-heave  model.  The  devices  also 
made  it  possible  to  conduct  two  special  tests  used 
to  validate  the  thaw  settlement  algorithm  of  the 
•  rost  heave  model 

Mathematical  model  of  frosl  heave 
and  thaw  settlement 

\  model  developed  under  an  earlier  cooperative 
lese.r.h  project  was  improved,  refined  and  ex¬ 
tended  I  he  refinements  were  introduced  as  verifi- 
.ation  work  progressed;  thev  included  accurate 
simulations  ol  pore  presstires,  temperatures  and 
t'ost  heave  \|j|or  extensions  to  (he  model  includ¬ 
ed  the  addition  ot  a  thaw  settlement  algorithm 
and  i  p-oktbilisiK  component  that  accounted  For 


variations  in  the  hvdrauhc  conductivity  caused  by 
spatial  and  temporal  variations  in  soil  properties. 

1  he  model  was  tested  bv  running  computer  sim¬ 
ulations  of  the  test  sections  at  Winchendon,  Mas¬ 
sachusetts,  and  at  Albany  County  Airport,  New 
\ork  I  he  calculated  frost  depths,  Frost  heaves 
and  thaw  settlements  were  compared  with  Field 
observations  at  the  test  sections 

In  general  the  results  show  that  good  results  can 
he  obtained  From  the  model  lor  soils  ranging  From 
silts  ti>  relatively  coarse-grained  and  only  margin 
ally  Frost  susceptible  soils  The  model  must  be 
calibrated,  however,  to  achieve  such  results.  This 
is  done  by  adjusting  the  value  For  hydraulic  con¬ 
ductivity  m  the  Freezing  zone  We  believe  that  the 
model  simulates  Freezing  processes  adequately  For 
engineering  purposes.  Model  outputs  include  Frost 
heave  at  the  pavement  surface  and  subsurface 
temperature  and  pore-water-pressure  distributions 
with  depth,  all  of  which  are  directly  useful  in  the 
pavement  design  process. 

Seasonal  variation  in  resilient 
modulus  of  granular  soils 

The  objective  of  this  phase  of  the  work  was  to 
develop  laboratory  methods  lor  characterizing  the 
seasonal  changes  in  the  resilient  modulus  of  gran¬ 
ular  soils  throughout  a  complete  annual  cycle  of 
Freezing,  thawing  (and  accompanying  loss  of  sup¬ 
porting  capacity  )  and  recovery  from  the  weakened 
condition.  A  laboratory  testing  method  was  devel¬ 
oped  For  conducting  repeated-load  triaxial  tests  on 
a  sample  initially  in  the  frozen  state,  continuing 
the  load  repetitions  after  the  sample  thawed,  and 
resuming  load  repetitions  again  at  successively  in¬ 
creased  levels  of  moisture  tension  obtained  by  pro¬ 
gressive  desaturation  of  the  sample. 

Tests  using  this  method  were  completed  on  1 1 
soils  From  paved  areas  at  test  sites  at  Winchendon, 
Massachusetts,  and  Albany  County  Airport,  New 
York.  By  means  of  multiple  regression  analysis  of 
the  test  results,  the  resilient  modulus  of  each  soil 
was  expressed  in  terms  of  the  governing  param¬ 
eters.  The  main  parameter  For  soils  in  the  frozen 
slate  is  temperature,  while  the  resilient  modulus  of 
soils  in  the  thawed  and  recovering  states  depends 
mainly  on  the  applied  stress  and  the  moisture  ten¬ 
sion. 

The  validity  of  the  method  was  tested  by  using 
the  expressions  For  resilient  modulus  to  calculate 
resilient  deflections  under  plate  loads  at  the  two 
test  sites  and  then  measuring  the  actual  deflections 
produced  by  a  Falling-weight  deflcclometcr  or  in 
some  cases  by  a  repeated-load  plate  hearing  de- 


'•icc.  Agreement  was  good  between  the  calculated 
and  measured  deflections  Accordingly,  it  is  con¬ 
cluded  that  the  laboratory  test  method,  supple¬ 
mented  bv  sotl-moisture-tension  and  teinpei attire 
data  obtained  from  field  measurements  or  predict¬ 
ed  by  the  trust -heave  model,  can  provide  evalua¬ 
tions  of  seasonal  variations  of  the  resilient  modu¬ 
lus  ot  granular  soils  and  base  courses. 

IMPLEMENTATION  OK 
RESEARCH  FINDINGS 

Investigations  in  the  principal  study  areas  have 
produced  results  that  can  be  advantageously  used 
in  designing  and  evaluating  pavements  in  frost 
areas.  The  examination  of  frost-susceptibility  in¬ 
dex  tests  (Chamberlain,  in  prep,  a)  identified  the 
Corps  of  Engineers  frost  design  soil  classification 
system  and  a  new  laboratory  freeze-thaw  test  as 
two  levels  of  testing  that  should  be  put  in  practice. 
The  frost -heave  model  can  be  implemented  benefi¬ 
cially  in  any  system  for  pavement  design  or  evalu¬ 
ation.  And  finally,  the  laboratory  repeated-load 
triaxial  test  on  thawed  and  recovering  soil  can  not 
only  be  implemented  in  mechanistic  design  or 
evaluation  systems,  but  when  used  with  either  the 
frost-heave  model  or  with  in-situ  measurements  of 
values  of  moisture  tension,  it  can  be  implemented 
in  systems  employing  a  cumulative  damage  ap¬ 
proach.  The  scope  and  extent  of  the  implementa¬ 
tion  of  each  of  the  research  findings,  and  its  po¬ 
tential  impact  on  pavement  design  and  evaluation, 
depend  on  the  type  of  system  that  is  in  use  for 
pavement  analysis. 

Corps  of  Engineers  frost  design 
soil  classification  system 

This  frost-susceptibility  classification  system 
(Table  1),  based  on  simple  classification  tests, 
leads  to  the  assignment  of  a  frost  group  number  to 
each  soil.  Its  implementation  (Fig.  41)  is  limited  to 
those  pavement  design  or  evaluation  systems  that 
are  based  on  the  frost  group  number. 

New  laboratory  freeze-thaw  test 

This  improved  freezing  test  incorporates  a  CBR 
test  on  the  thawed  specimen.  Outputs  (Fig.  42)  are 
heave  rate  from  the  freezing  phase,  and  CBR  in 
the  thawed  state;  scales  of  susceptibility  tn  frost 
heave  and  to  thaw  weakening  may  be  derived  from 
these  outputs  (Table  3).  Frost-heave  susceptibility 
can  be  applied  effectively,  as  an  indication  ot  po¬ 
tential  roughness,  to  any  mechanistic  or  empirical 


design  or  evaluation  system.  The  classification  of 
susceptibility  to  thaw  weakening  can  be  used  as  an 
indicator  of  seasonal  (spring)  loss  of  support  as 
part  of  any  empirical  design  or  evaluation  system. 
In  developing  a  pavement  design  these  indicators 
would  be  used  as  adjustments  to  design  thickness¬ 
es  determined  by  other  types  of  analyses. 

Frost-heave  model 

The  frost-heave  model  provides  for  the  first 
time  an  ability  to  calculate  with  reasonable  confi¬ 
dence  the  magnitude  of  the  heave  that  can  be  ex¬ 
pected  in  a  given  pavement  cross  section  under 
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Figure  43.  Implementation  of  the  frost-heave  model. 


prescribed  climatic,  geotechnical  and  hydraulic 
conditions.  Since  the  outputs  of  depth  of  frost  and 
magnitude  of  heave  (Fig.  43)  are  referenced  to  a 
particular  point  on  the  pavement  where  conditions 
are  known,  the  model  does  not  predict  pavement 
roughness.  Heave  at  a  point  can  be  used  as  an  in¬ 
dicator  of  potential  roughness,  however,  and  can 
be  implemented  as  an  adjunct  to  any  pavement  de¬ 
sign  or  evaluation  system.  For  example,  the  calcu¬ 
lated  frost  heave  might  serve  as  a  basis  for  adjust¬ 
ing  a  trial  design  thickness  if  necessary  to  reduce 
the  expected  winter  pavement  roughness.  The  sec¬ 
ond  principal  output  from  the  model,  the  predict¬ 
ed  depth  of  frost  beneath  a  pavement  having  a  cer¬ 
tain  trial  cross  section,  can  be  used  as  direct  input 
for  those  design  systems  in  which  the  total  thick¬ 
ness  of  the  pavement  section  depends  on  the  depth 
of  frost.  It  can  also  be  used  as  an  adjunct  to  any 
design  system. 


mented  directly  in  any  mechanistic  design  or  eval¬ 
uation  system  that  employs  a  multilayered  or 
finite-element  simulation  model  formulated  to 
analyze  nonlinear  materials. 
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Repeated-load  triaxial  test 
on  frozen  and  thawed  soil 

The  repeated-load  triaxial  test  on  frozen  and 
thawed  specimens  provides  a  means  of  evaluating 
the  resilient  modulus  of  subgrade  and  base  soils  at 
various  stages  during  the  freeze-thaw-recovery 
cycle.  The  regression  equations  for  soil  in  the 
thawed  state  (Fig.  44)  are  of  the  greatest  interest, 
as  in  many  cases  they  represent  the  condition  hav¬ 
ing  the  lowest  resilient  modulus  and  consequently 
the  greatest  potential  for  pavement  distress.  The 
expressions  for  resilient  modulus  can  be  imple- 
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Figure  44.  Implementation  of  the  repeated-load  tri¬ 
axial  test  for  stress-strain-deflection  analysis. 
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Figure  45.  Implementation  of  the  seasonally  varying  resilient 
modulus  for  stress-strain-def lection  analysis  with  a  cumulative 
damage  approach. 


Evaluation  of  seasonal  variation 
of  resilient  modulus 

It  is  unrealistic  to  base  the  design  or  evaluation 
of  a  pavement  on  only  the  lowest  resilient  modulus 
value  reached  during  the  year  (usually  during 
thawing),  and  there  is  no  reasonable  basis  for  se¬ 
lecting  any  other  single  value  to  serve  as  an  annual 
average  representative  of  all  the  seasons.  Rather, 
methods  that  include  a  cumulative  damage  ap¬ 
proach.  currently  coming  into  greater  use,  offer 
the  advantage  of  more  rigorously  assessing  the  ef¬ 
fect  of  the  complete  annual  cycle  of  freezing, 
thawing  and  recovery.  Application  of  these  meth¬ 
ods  requires  that  the  resilient  modulus  be  ex¬ 
pressed  as  a  function  of  time,  so  that  the  pavement 
performance  can  be  analyzed  by  dividing  the  year 
into  discrete  intervals  during  which  the  modulus 
may  be  assumed  to  be  constant.  The  evaluation  of 
seasonal  variations  of  the  resilient  modulus  re¬ 
quires  that  the  modulus  be  characterized  in  terms 
of  the  moisture  tension  and  that  the  seasonal  vari¬ 


ation  of  the  moisture  tension  be  monitored  or  pre¬ 
dicted.  Figure  45  shows  that  repeated-load  triaxial 
tests  can  be  used  to  measure  the  resilient  modulus 
as  a  function  of  temperature,  soil  stress  and  mois¬ 
ture  tension.  The  frost-heave  model  and  its  associ¬ 
ated  thaw-settlement  model  are  used  to  predict  the 
key  parameters  of  temperature  and  moisture  ten¬ 
sion  as  variables  in  time  and  space  With  this  link 
age  the  resilient  modulus  at  various  depths  is  de¬ 
fined  as  a  continuous  function  of  time,  facilitating 
the  application  of  mechanistic  analyses  lor  pa\e 
ment  design  and  evaluation  using  a  cumulative 
damage  approach. 
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